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Chapter 1 
Introduction 



The Standard Model of particle physics is a successful theory describing the elementary 
particles and their interactions. It combines the 811(2)^(8)11(1)^ theory of electroweak 
interactions with the SU(3) Quantum Chromodynamics (QCD) - the theory of strong 
interactions. The strong interacting particles, hadrons (e.g. protons), are build from 
more elementary point-like objects, quarks and gluons. The quarks are charged and 
they interact also electromagnetically with other charged particles, e.g. electrons, by 
exchange of the point-like photons. The electromagnetic interactions are described 
by the Quantum Electrodynamic (QED), being a part of the unified SU(2)i(8)U(l)y 
theory. 

This work is devoted to the theoretical study of the photoproduction of the photon 
with a large transverse momentum, called the Deep Inelastic Compton (DIG) process, 
in the electron-proton scatterings at the DESY HERA collider. We calculate the 
cross sections for this process with the beyond leading logarithmic accuracy within a 
framework of the perturbative QCD (see eg. [1]). 

The photoproduction is a process in which a real or almost real photon collides 
with another particle producing some final state particles (for recent reviews see [2, 
3, 4]). Currently, all high energy photon-proton collisions are realized in the electron- 
proton scatterings. They correspond predominantly to events with low virtuality of the 
exchanged photons, ^ 0. These quasi-real photons colliding with the proton may 
lead for instance to a production of jets or particles with a large transverse momentum, 
Pt- If Pt is much larger than the QCD scale, px ^ ^qcd, then the perturbative QCD 
can be applied to describe such a process. 

Processes in which the large transverse momentum of jets, hadrons or photons is 
observed in a final state in the e"'"e~, hadron-hadron and electron-proton colliders play 
an important role in testing QCD and in measuring the parton densities or the parton 
fragmentation functions as well as the strong coupling constant as- The process with 
the production of the photon is a special one, since the photon may couple directly 
to quarks involved in the hard QCD process, so it may provide a relatively clear 
information about the QCD dynamics. Such photons arising predominantly from the 
hard process are called prompt photons. Although the cross sections for the prompt 
photon production are smaller than for the production of jets, these processes are 
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considered as an important source of complementary informations. 

The photoproduction of the prompt photon in the electron-proton scattering is 
shown in Fig. 1.1. In such a reaction an almost real photon emitted from the electron 




Figure 1.1: The ep ^ ej X photoproduction. 



interacts with the proton leading to the production of large-pr photon and a hadronic 
system X: ep — *■ ejX ^. Particles in the hadronic final state can form a jet or jets 
which balance the photon transverse momentum. If, beside the final photon, the jet is 
considered, we write: ep — 67 jet X (Fig. 1.2), where X stands for other hadronic final 
state. In the work we study both types of processes: ep — 67 X and ep — >■ 67 jet X. 




Figure 1.2: The ep —>■ 67 jet X photoproduction. 

In the lowest order of perturbative QCD the photon emitted by the electron is 
scattered from the quark being a component of the proton (the Born contribution), 
see Fig. 1.3. The final state of the hard interaction between the photon and the initial 
quark consists of the large-p-p photon and the quark from which the jet arises. In this 
lowest order process both the mediating and the final photon interact directly with 
the quark. 

""^For a relatively small momentum transfer between the initial and the final electron, « 0, the 
exchange of the photon dominates over the Z boson exchange and we do not include the Z boson in 
our analysis. We also neglect the emission of the final photon directly from the electron as it gives a 
very small contribution to the cross section for the photoproduction of large-p^ photons [5] . 
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Figure 1.3: The Born contributions to the ep —>■ e^{jet) X photoproduction. 



However, there are also processes of a quite different nature. At high energies the 
photon being a point-hke gauge boson, may interact hke a hadron. This hadron-hke 
properties of the photon are due to fluctuations into a virtual pair quark-antiquark 
or into a vector meson, see Fig. 1.4. Due to these fluctuations the photon may 




Figure 1.4: Various processes contributing to the hadron-hke "structure of the photon". 



behave as an object consisting of quarks, antiquarks and gluons. The interactions 
in which such a photon, called the resolved photon, exhibits its complex hadronic- 
like "structure" are effectively described using the formalism of the photon structure 
functions [6, 7] in full analogy to the corresponding formalism for the proton [8, 9, 1]. 
The density of probability of "finding" the parton (quark or gluon) "inside" the photon 
is given in this formalism by the parton density/distribution in the photon, fq{g)/-y 
The process 'J —>■ q + anything described by the function fg/^ is illustrated in Fig. 
1.5. The parton densities are related to the structure functions, as for the proton. 



^AAAA/C^J■ 



Figure 1.5: The resolved photon. 

However, the photon structure functions, e.g. F2 , in contrast to the corresponding 
proton structure functions, are calculable in the Parton Model, based on the 7*7 — > qq 
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process. The F2 is proportional to the electromagnetic coupling constant, a, and it 
depends logarithmically on the energy scale of the process, /x, already in the Parton 
Model 2; 

F^^a\n-^ + ... (1.1) 

We stress that this logarithmic scaling violation in F2 in the Parton Model arises 
from the purely electromagnetic coupling 7 — > qq. For more detailed discussion on the 
concept of the "photon structure" see e.g. [3, 4, 10]. 

So, taking into account the hadronic structure of the photon, the photon mediat- 
ing in the electron-proton scattering shown in Figs. 1.1 and 1.2 may interact with the 
constituents of the proton directly (as in Fig. 1.3) or indirectly as the resolved photon. 
Similarly, the final photon can be produced directly in the hard interaction with the 
quark (as in Fig. 1.3) or it may originate from fragmentation processes. The cross 
section of the latter processes involve relatively poorly known parton-to-photon frag- 
mentation functions (see e.g. [11] and references therein) which describe the density of 
probability of the parton "decay" into the photon, see Fig. 1.6. Fortunately, one can 



Figure 1.6: The quark-to-photon fragmentation. 

suppress this contribution introducing the isolation of the final photon. The isolation 
constraints, which separate the photon from energetic hadrons, are introduced in both 
experimental and theoretical analyses. We follow this approach in this work. 

The standard experiments to extract the parton densities in the proton or pho- 
ton are the Deep Inelastic Scatterings in the ep or cy collisions, where the structure 
functions and F2 are measured. From these structure functions the correspond- 
ing parton densities can be derived. The photoproduction processes, as the prompt 
photon production considered herein, can provide an additional information about the 
parton densities in both the proton and the photon, see e.g. [2, 3]. 

The production of the photons with large transverse momenta in the 7p scatterings, 
and its importance for probing the proton structure, was studied theoretically for 
the first time in 1969 [8]. The concept of the resolved photon was included in the 
analyses of this process ten years latter [12]. Then, a few groups of authors studied 
the inclusive (i.e. non-isolated) large-pT photon production in the jp or ep scattering, 
and its importance for testing QCD and constraining the parton densities in both 
the proton and photon [13]- [26]. The next-to- leading order (NLO) QCD predictions 
for the photoproduction of isolated photons in the ep process at the HERA collider 

^In the case of the 7*7 — > qq process, the energy scale /i squared is equal to the virtuality, Q^, of 
the virtual photon (7*) which probes the structure of the real photon (7). 
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have been presented in the papers [27, 28, 29, 30]. Finally, the photoproduction of 
the isolated photon and a jet at HERA has been calculated in NLO QCD in [31]- 
[36] . There is also a NLO calculation for the isolated photon and the isolated photon 
plus a jet production in the deep inelastic events (Q^ 3> 1 GeV^) [37, 38, 11, 39], 
which is in many theoretical aspects close to the calculations for the corresponding 
photoproduction events. 

The photoproduction of isolated photons without and with additional jets has been 
measured at the HERA collider by the ZEUS [40]-[45] and HI [46]-[49] Collaborations. 
All the NLO QCD predictions as well as the Monte Carlo simulations for the isolated 
photon with no jet requirement tend to lie below the ZEUS data for the photon ra- 
pidities 1]"' < 0.1 [41, 42, 43, 45], and below the final HI data [49] in the whole range 
of the photon rapidities On the other hand the HI data [49] for the isolated photon 
plus jet production are somewhat better described by the QCD predictions. 

In the paper [44] the ZEUS Collaboration has implemented in Monte Carlo simu- 
lations for the prompt photon plus jet production the intrinsic transverse momentum 
of partons in the proton. From a fit to the data it was found that the effective intrinsic 
transverse momentum is very large with respect to the mass of the proton, namely 
< kr >— 1.69 ± 0.18^0 20 GeV. In our work the intrinsic transverse momentum is 
not included, since one can argue that the measured < /c^ > in the proton describes 
effectively higher order emissions or multiple interactions between particles involved 
in the process. Moreover, it was shown in [32] that the NLO QCD calculation is able 
to describe the ZEUS data [44] with no need for the intrinsic transverse momentum. 

As it was mentioned above, there are some moderate discrepancies between the 
QCD predictions and the data for the photoproduction of the isolated photons in ep 
scattering, especially when no jet in the final state is considered. Note, that some 
differences between theory and data are also observed for other processes involving 
photons in the hadron-hadron, ep and e+e~ collisions, see e.g. [50, 51, 52, 53]. It may 
indicate, that our understanding of such processes is not satisfactory yet, and further 
theoretical and experimental searches are needed. This is one of reasons for our study. 

This thesis contains a detailed theoretical analysis of the cross section for the 
photoproduction of the prompt photon and the prompt photon plus a jet including 
effects due to the structure of the photon and the fragmentation into photon. In 
most of theoretical studies the parton densities in the photon, and the parton- 
to-photon fragmentation functions, D^, are treated as quantities of order 0{a/as). 
In our calculations they are considered as quantities of order 0{a). This leads to a 
different set of diagrams included in our calculation [28, 30, 33] and in calculations 
of other authors [27, 29, 31, 32, 34, 35, 36] for the photoproduction of the isolated 
photons at HERA. 

We implement the isolation restrictions in the cross section using two methods. 
First, we use the small cone approximation method, in which the isolation is imple- 
mented in an approximated way. This method was previously used by other authors. 

^The rapidity of e.g. the photon is defined as i)'^ = — lntan(0/2), where 9 stands for the angle 
between the final photon momentum and the momentum of the initial proton. The positive rapidity 
is pointed in the proton direction. 
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However, not all expressions needed to calculate the higher order corrections in the 
small cone approximation exist in the literature, so we derive and present the missing 
formulae. Next, we use the method, in which the phase space of produced particles 
is divided into a few parts. This allows to implement the isolation, as well as other 
kinematic cuts, in an exact way. We perform a comparison between both methods 
and show that the small cone approximation is quite accurate and leads to reliable 
predictions for the prompt photon production at HERA. 

The division (slicing) of the phase space is a standard approach to calculate various 
cross sections, e.g. cross sections for the isolated photon production. However, our 
method differs in details from methods applied in other calculations and allows to 
obtain relatively simple analytical formulae for the cross sections including higher 
order corrections. 

The formulae for the higher order corrections are different in each part of the 
divided phase space. Our method of the division allows to apply in some parts of the 
phase space the known cross sections for the non-isolated photon production. In other 
parts of the phase space we derive all needed analytical expressions for the higher 
order corrections to the isolated photon production in the ep collision. Some obtained 
expressions are consistent with formulae presented previously by other authors (for 
processes other than this considered herein). We find it useful to present in this work 
in a compact analytical form all the expressions, which are necessary to obtain the 
cross section for the isolated photon production in the ep scattering including higher 
order corrections. 

This thesis is based on the papers listed below and it contains in addition some 
new results never published before: 

• M. Krawczyk and A. Zembrzuski, "Probing the structure of virtual photon in 
the deep inelastic Compton process at HERA," Phys. Rev. D 57 (1998) 10 
[arXiv:hep-ph/9708274], 

• M. Krawczyk, A. Zembrzuski and M. Staszel, "Survey of recent data on photon 
structure functions and resolved photon processes," DESY 98-013, arXiv:hep- 
ph/9806291, 

• M. Krawczyk and A. Zembrzuski, in: A. Astbury, D. Axen, J. Robinson (Eds.), 
Proceedings of the 29th Int. Conference on High Energy Physics, ICHEP'98, 
Vancouver, Canada, July 1998, World Scientific, 1999, p. 895, "NLO prediction 
for the photoproduction of the isolated photon at HERA," arXiv:hep-ph/9810253, 

• M. Krawczyk and A. Zembrzuski, "The forward photon production and the 
gluonic content of the real and virtual photon at the HERA collider," Nucl. 
Phys. Proc. Suppl. 82 (2000) 167 [arXiv:hep-ph/9912368], 

• M. Krawczyk, A. Zembrzuski and M. Staszel, "Survey of present data on photon 
structure functions and resolved photon processes," Phys. Rept. 345 (2001) 265 
[arXiv:hep-ph/0011083]. 
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• M. Krawczyk and A. Zembrzuski, "Photoproduction of the isolated photon at 
DESY HERA in next-to-leading order QCD," Phys. Rev. D 64 (2001) 114017 
[arXiv:hep-ph/0105166], 

• A. Zembrzuski and M. Krawczyk, "Photoproduction of isolated photon and jet 
at the DESY HERA," arXiv:hcp-ph/0309308, 

The paper is organized as follows. First, in Chapter 2 we discuss our choice of 
diagrams and present some calculation details. The formula for the cross section and 
our results for the inclusive (non-isolated) photon production are briefly discussed in 
Chapter 3. In Chapter 4 the isolation cuts are defined and the calculation for the 
isolated photon in the small cone approximation is presented; we study the influence 
of the isolation cut on the production rate of the photon and the role of other cuts 
applied in experiments. Then, in Chapter 5 the division (slicing) of the three body 
phase space is discussed and we present our analytical results. Here, the numerical 
predictions for the isolated photon production with the exact implementation of the 
isolation cuts are also presented and compared with the previous approximated ones. 
The results for the ep 67 jet X process are presented in Chapter 6. We discuss 
theoretical uncertainties and compare our predictions with predictions of other NLO 
QCD calculations, as well as with existing data for the photoproduction of the isolated 
photon or the isolated photon plus a jet at HERA (Chapters 4, 5, 6). In Chapter 7 
we briefly discuss the sensitivity of the leading order (LO) cross section to the gluon 
distribution in the photon taking into account a non-zero virtuality of the exchanged 
photon, 7^ 0. The summary is given in Chapter 8. Finally, the Appendices contain 
the formulae for the cross sections for all the processes included in our analysis. Some 
of these formulae are taken from the literature and the other are briefly derived. 
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Chapter 2 



The Deep 



Inelastic Compton 
NLO QCD 



process m 



Our calculation bases on the standard leading twist perturbative QCD description 
of hard hadronic processes. This means that we consider such diagrams where one 
active parton from each initial particle is involved in the hard process. However, in 
some Monte Carlo simulations performed by experimental groups recently (2004) a 
large effect of multiple interactions was found for the prompt photon production at 
HERA. This effect is included in experimental analyses by introducing some models. 
We present predictions based on our calculations and "corrected" by an experimental 
group for the multiple interactions as well as hadronization processes in Chapters 5 
and 6. 

There are several calculations for the photoproduction of prompt (isolated) pho- 
tons at the HERA collider, namely: the calculation of Gordon and Vogelsang (GV) 
[27], Gordon (LG) [31], Krawczyk and Zembrzuski (K&Z) [28, 30, 33] and Fontannaz, 
Guillet and Heinrich (FGH) [29, 32]. All these calculations differ from one another by 
set of diagrams included in next-to-leading order (NLO) . Below we discuss our choice 
of diagrams (Sec. 2.1), as well as some calculation details (Sees. 2.2, 2.3). 

2.1 Contributing processes 

2.1.1 Born process and 0{as) corrections, resolved 7 or frag- 
mentation into 7 

We wish to study the hard electron-proton scattering leading to a production of a 
photon or a photon plus a jet and anything else, ep 67 {jet) X (Figs. 1.1-1.3). In 
this reaction the mediating photon arising from the electron interacts with a partonic 
constituent (quark or gluon) inside the proton [8]. In the lowest order (Born) process, 
7? ~^ IQ1 fhe photon is scattered from a quark yielding the final photon and the quark, 
see Fig. 2.1. This subprocess is of order 0{a^) and has a pure electromagnetic nature. 
To obtain the NLO QCD predictions, the corrections of order 0{as) to the Born 
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process have to be taken into account. These corrections lead to partonic processes 
of order 0{a'^as). They include the virtual gluon exchange (Fig. 2.2), the real gluon 
emission (Fig. 2.3) and the process •yg — > •yqq (Fig. 2.4) [15, 17, 18, 25]. 

The corrections due to the virtual gluon exchange and the real gluon emission con- 
tain infrared singularities which cancel when both contributions are added up properly 




+ crossed diagram 



Figure 2.1: The Born process (the Compton scattering on the quarlc). 




Figure 2.2: An example of the virtual gluon corrections to the Born process. 
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Figure 2.3: An example of the real gluon corrections to the Born process, jq 'jqg- 




Figure 2.4: An example diagram for the process yg — > yqq. 
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in the scattering amplitude squared (see Sees. 3.1, 4.2, 5.2). The contribution due 
to the processes •yq •yqg and jg ygq contains another type of singularities, so 
called mass or collinear singularities, which do not cancel. In order to remove them 
from the cross section the factorization procedure is applied: the singularities are sub- 
tracted and shifted into corresponding parton densities in the proton or photon or 
into parton-to-photon fragmentation functions [18, 25]. At this stage the factorization 
scale appears. The bare {scale invariant) parton densities in the proton are replaced 
by renormalized scale dependent densities. Moreover, the scale dependent parton den- 
sities in the photon and the scale dependent fragmentation functions appear in the 
calculation. They are necessary ingredients in the NLO calculation to absorb the mass 
singularities. 

The factorization procedure is illustrated in Figs. 2.5 and 2.6: the singularities 
due to the collinear configurations in the vertexes j ^ qq and q ^ ■jq are shifted and 



collinear singularity 




Figure 2.6: A sample of diagrams illustrating the factorization of the collinear singularities 
from the process "yq — > ^qg into the fragmentation function, D^/q- 
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absorbed by the corresponding quark density in the photon (/ij/-y) or by the quark-to- 
photon fragmentation function (-D^/g) , respectively (for precise definitions and details 
see [18, 25]). 

It is worth mentioning that there was a discussion whether for the photoproduc- 
tion of isolated photons in e^e~ collisions the conventional factorization breaks down, 
and whether the cross section is an infrared safe quantity [54, 55, 56]. In princi- 
ple these questions could as well occur for the photoproduction of isolated photons 
in ep collisions. However we do not deal with this problem because it arises from 
2 — i> 3 processes involving the parton-to-photon fragmentation, which are absent in 
our calculation. We checked this explicitly and found that all the singularities in our 
calculations for the isolated photon production are canceled or factorized, as in the case 
of non-isolated photon production, and the cross sections are well defined [28, 30, 33] 
(see also [27, 29, 31, 32], [34]-[39]). 

The processes involving the parton densities in the resolved initial photon are 
shown in Fig. 2.7. In these processes a parton (quark or gluon) from the photon 
interacts with a parton arising from the proton yielding the direct final photon and a 
parton. Fig. 2.8 shows the processes with the fragmentation of a parton into the final 
photon. Here the direct initial photon interacts with a parton from the proton leading 
to a production of two partons; one of these partons produces the final photon in the 
fragmentation process. The name direct photons stand for the photons participating 
directly in the hard partonic process. 




Figure 2.7: The processes with the resolved initial photon. 
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Figure 2.8: The processes with the parton-to-photon fragmentation. 

2.1.2 Counting of orders. Alternative set of diagrams 

The hard partonic processes in Figs. 2.7 and 2.8 are of order 0{aas)- They are 
convoluted with the corresponding parton densities or fragmentation functions which 
are proportional to the electromagnetic coupling constant a. These distributions con- 
tain the logarithmic dependence on the factorization/renormalization scale /i: ~ 
aln/i^/Ag^^ and ~ aln/i^/Ag^^. The logarithm behaves like the inverse of the 
strong coupling, In/i^/Ag^^ ~ l/a5 and many authors (e.g. [27, 29, 31, 32, 34, 35, 36]) 
treat and as being of order 0{a/as)- This conclusion arises from the evolution 
equation and the renormalization group equation, as it shown in [57]. If so, then the 
processes shown in Figs. 2.7, 2.8 are of order 0{aas-a/as)= 0{a'^), i.e. identical with 
the order of the Born process (Fig. 2.1). Note however, that the logarithmic behavior 
of the parton densities in the photon or the parton-to-photon fragmentation functions 
arises from the pure electromagnetic coupling •y ^ qq and q •yq, respectively [6], 
as already discussed in Introduction. From this point of view the considered distri- 
butions should be rather treated as the quantities of order 0{a) since the logarithm 
ln/i^/AQ(-7£) coming from the electromagnetic interaction is not related to the strong 
coupling constant, see also [20, 28, 30] and for more detailed discussion in [58]. Such 
a counting, and ~ 0{a), leads to an assignment of the 0{a'^as) order to the 
processes shown in Figs. 2.7, 2.8. 

The different counting of orders in the strong coupling leads to the different sets of 
processes included in our NLO calculation for the prompt photon production at HERA 
[28, 30, 33] in comparison with other NLO calculations. They include in addition 
(beside the 0{as) corrections to the Born cross section) the 0(0:5) corrections to 
the processes with the resolved initial photon or parton-to-photon fragmentation [27, 
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29, 31, 32], and the 0{as) corrections to the processes with both the resolved initial 
photon and parton-to-photon fragmentation [27, 29, 32] ^. These corrections are not 
included in our NLO calculation since in our approach they should be taken into 
account together with the next-to-next-to-leading (NNLO) order corrections to the 
Born process. 

It should be emphasized that the 0{as) corrections to the Born process contain 
(after the factorization procedure) an explicit logarithmic dependence on the scale /x, 
Ins///^. This dependence is compensated by In h'^/Aq^jj from and if the 0{as) 
corrections to the Born process are included in the cross section consistently with 
the corresponding contributions involving parton densities in the photon or parton- 
to-photon fragmentation functions. Similarly, the 0{as) corrections to the processes 
with the resolved initial photon or the fragmentation into the final photon should be, 
from our point of view, taken into account together with the 0{a^) corrections to the 
Born process in order to compensate the dependence on the /i^ in the hadronic cross 
section. 

One can argue that, since the other authors [27, 29, 31, 32] include more diagrams, 
their calculation is more accurate. It is true that including more diagrams usually 
improves the quality of QCD predictions but it is not obvious that it is always the 
case (see also the discussion in Sees. 4.4.3 and 4.4.5). 

It is worth mentioning that there are also other authors who treat the processes 
involving photons in a similar way as we do. The authors of [64] (sec also [65, 66]) 
present the NLO calculation for the prompt photon production in the e+e^ collision 
with the quark-to-photon fragmentation. They include the fragmentation function at 
the same order as the 0{as) corrections, since the fragmentation is in a close relation 
with the corresponding coUinear configuration in 0{as) corrections (compare Fig. 2.6). 
It allows to cancel an explicit dependence on the scale /x^. A similar NLO calculation 
for the prompt photon production in the deep inelastic ep events (with > 10 GeV^) 
at the HERA collider is presented in [39, 38, 11]. In this calculation the diagrams shown 
in Figs. 2.1-2.4 and 2.8 are included without the 0{as) corrections to the processes 
involving the fragmentation (the contributions due to the resolved initial photon are 
negligible for large Q^). The cancellation (to a large extent) of the //^ dependence is 
therein also stressed. 

Another kind of arguments is presented in [67] , where the NLO calculation for the 
prompt photon production in the e'^e~ collision is investigated. The authors of [67] 
claim that the parton-to-photon fragmentation, despite being of order 0{a/as) for 
the non-isolated photon production, should be counted as the quantity of order 0{a) 
for the isolated final photon, since the isolation itself is a correction of order 0{as)- 

Counting and as being of order 0{a/as) allows for a self-consistent expansion 
of physical quantities (structure functions, cross sections) in powers of I/Iu/x^/Aq^^ 
(with as ~ 1/ Iu/x^/Aqp^) -|- ...) [57]. On the other hand, our approach, where f^, 
Dy ~C(q;), allows for a consistent counting of powers of as (since the logarithm 
Iuii^/Aqq^ in /-y and it is not 1/as, see also [20, 58]). 

^AU these QCD corrections have been calculated (most of them twice) in [59]- [63] 
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The photon is a very special particle: being the point-like object it sometimes 
exhibits hadronic-hke "structure". This double nature of the photon leads to more 
complicated description within the QCD than for purely hadronic processes, and fur- 
ther study is needed to clarify what is the proper organization of the QCD perturbative 
series for processes involving photons. This should include among others also the cal- 
culation of the 0{al) corrections to the Born process. 

2.1.3 Full set of diagrams included in analysis. Resolved 7 
and fragmentation into 7, box diagram 

Besides the diagrams discussed in Sec. 2.1.1, we include in our numerical analysis also 
other diagrams, namely the diagrams with the resolved 7 and the fragmentation into 7 
and the box diagram. In Figs. 2.9, 2.10 the processes with the resolved initial photon 
and the fragmentation into the final photon are shown. If we take ~0(q;), 
all these processes are of order 0{a'^a'g) and strictly speaking go beyond the NLO 
accuracy of our calculations. Nevertheless we include them, since they were taken 
into account in most of existing NLO calculations for the (non-isolated or isolated) 
prompt photon production at HERA and are found to be important [15, 18], [20]-[24], 
[26]-[34]. 

Finally, we take into account the photon-by-gluon scattering, so called box process, 
19 ~^ 19^ shown in Fig. 2.11 [68]. The box process is also of order (9(a^Q;|), i.e. beyond 
the accuracy of our calculation. However it is a very special process as it is the lowest 
order contribution to the photon-by-gluon scattering, it contains no singularities, and 
no parton densities in the photon or fragmentation functions are involved. We think 
that including the box process does not introduce any additional ambiguities in the 
summation of partonic cross sections. The box contribution to the reaction ep e'-fX 
or 7p — i> 7X was considered in previous calculations and it is known to be large [14, 18], 
[19]-[24], [26]. It is included in the calculations for the isolated photon production at 
HERA by authors K&Z [28, 30, 33] and FGH [29, 32] but it is omitted by GV [27] 
and LG [31]. 

To summarize, in our NLO calculation we take into account the following contri- 
butions: 

• the Born contribution (Fig. 2.1), 

• the 0{as) corrections to the Born diagram, including the jg — > ■yqq process 
(Figs. 2.2-2.4), 

• the processes with the resolved initial photon (Fig. 2.7) or with the fragmentation 
into the final photon (Fig. 2.8), 

• the processes with the resolved initial photon and the fragmentation into the 
final photon (Figs. 2.9, 2.10), 

• the box process (Fig. 2.11). 
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Figure 2.10: Continuation of Fig. 2.9. 




+ crossed diagrams 



Figure 2.11: The box diagram. 



2.2 Calculation details 

We perform NLO QCD calculations in the modified Minimal Subtraction (MS) renor- 
malization scheme [69]. The factorization/renormalization scales in parton densities 
and fragmentation functions are assumed being equal to the renormalization scale in 
the strong coupling constant and are denoted as /i. As a reference we take fi equal 
to the transverse momentum ^ (transverse energy) of the final photon, ^ = Pr = Ej,. 

^By the transverse momentum, p^, we mean the component of the momentum perpendicular to 
momenta of the initial particles. 
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For a comparison = Ej,/2 and n = 2Ej. will be also considered. 

The quark masses are neglected in the calculation and the number of active flavors 
is assumed to be Nf—A or, for a comparison, Nf—3 or 5. The cross section is propor- 
tional to the electric charge of the quark in the fourth power, so the contribution of 
the bottom quark (e;,=-l/3) is expected to be much smaller than the contributions of 
the up and charm quarks (e„=ec=2/3), and the predictions for iVj=4 and 5 should not 
differ considerably. On the other hand, the differences between the results obtained 
using Nf—A and 3 can be large. 

The two-loop coupling constant ag is applied in the form 

^ .0 ^ 127V _ 6(153 - hi[hi(/,7A|,,J] 

'^^^ {33-2Nf)ln{^ly^lcny (33 - 2Ar^)2 M/^VA^cd) ^' ^'^ 

To obtain the QCD parameter Aqcd appearing in the strong coupling constant (2.1), 
we use the world average oi as at the scale Mz (the mass of the Z° gauge boson) [70] : 

a^(M|) = 0.1183 ± 0.0027. (2.2) 

In order to minimize theoretical and experimental uncertainties the above as{Mz) 
value was determined in [70] from precise data based on NNLO analyses only; the 
data given at scales different than Mz were extrapolated to the Mz scale using the 
four-loop coupling. Although we use the two-loop expression (2.1), we apply (2.2) 
as the best estimation of the true value of as{Mz). We take the number of active 
flavours Nf=3, 4 or 5 at scales n < nic, rric < fJ. < nib aiid < ^^fc; respectively, with 
the following charm, bottom and Z masses: rric — 1.5 GeV, = 4.7 GeV, Mz — 91.2 
GcV [70]. Assuming that as{^'^) is a continuous function at /x = and fi = nib the 
obtained Aqcd parameters are: Aqcd=0-386, 0.332 and 0.230 GeV for Nf=3, 4 and 
5, respectively. The above Aqcd values are used in numerical calculations discussed 
in the next chapters (Chapters 3, 5, 6) ^. 

In the calculations we apply the Gliick-Reya-Vogt (GRV) parton densities in the 
proton [71] and photon [72], and the GRV fragmentation functions [73]. For a com- 
parison we also use other parametrizations, namely Martin-Roberts-Stirling-Thorne 
(MRST98) [74], (MRST99) [75], (MRST2002) [76], CTEQ4M [77], CTEQ6M [78], Au- 
renche-Chiappetta-Fontannaz-Guillet-Pilon (ACFGP) [24], Aurenche-Guillet-Fontan- 
naz (AFG) [79] and (AFG02) [80], Gordon-Storrow (GS) [81], Cornet- Jankowski-Kraw- 
czyk-Lorca (CJKL) [82], Duke-Owens (DO) [15], Bourhis-Fontannaz-Guillet (BFG) 
[83] and Gliick-Reya-Stratmann (GRS) [84]. 

Following experimental analysis [40]- [49], we consider the photoproduction of the 
photon at HERA with transverse momentum (transverse energy) higher than 5 GeV. 
In such processes the emission of the final large-p^ photon directly from the electron 
(Betlio-Hoitler ])i()( oss) is negligible [5] and we omit it ^ (this process is neglected in all 

^Note that in Chapters 4 and 7 different numbers are used for consistency with the published 
results. 

*For transverse energy higher than 5 GeV the momentum transfer, -t (see Eq. (A. 6)), is higher 
than 25 GeV^. Note that for Q^^x = 1 GeV^ and lower values of the momentum transfer, 4< —t < 10 
GeV^, the contribution of the Bethe-Heitler process may be no-negligible in some kinematic regions, 
see Fig. 6 in Ref. [85]. 
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existing calculations for the pliotoproduction of isolated photons at HERA [27]- [34]). 

The results presented in next sections are obtained in NLO QCD with use of the 
GRV set of parametrizations with = Ej., Nf = 4 and Aqc'd=0.332 GeV unless 
stated otherwise. 



2.3 Equivalent photon approximation 

Our aim is to consider the production of photons with large transverse momentum in 
the electron-proton scattering in processes in which the electron is scattered at a small 
angle. In such events the mediating photon is almost on-shell, ~ (photopro- 
duction events), and the cross section can be calculated using the equivalent photon 
(Weizsacker- Williams) approximation [86, 87] (see also e.g. [88, 89, 4]). In this approx- 
imation the differential cross section for the ep collision is related to the corresponding 
differential cross section for the 7p collision; in case of the considered herein reaction, 
ep — > ejX (or ep 67 jet X), we get (see Fig. 2.12): 

^^ep^e^Uet)x ^ J Q ^ / ,{y)da^P^^^^^'^^ dy , (2.3) 

where y is the fraction of the initial electron momentum carried by the exchanged 
photon, and G-^/e stands for the flux of the real photons emitted from the electron. 
We use the photon spectrum in the form [89]: 

Gi/e{y) = Tr\ n — — y^~y~iT2 — i' ^.4 

27r y mly^ y Q^„^. 

with irtf, being the electron mass. In the numerical calculations the maximal photon 
virtuality Ql^ax = 1 GeV^ is assumed, what is a typical value for the recent photopro- 
duction measurements at the HERA collider [40]- [49]. The above formula describes 
the spectrum of equivalent real (transversally polarized) photons. We do not take into 
account longitudinally polarized photons and the interference between longitudinally 
and transversally polarized photons, since they give a very small contribution in the 
kinematic regions which we study [5, 90]. 
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Chapter 3 

Inclusive photon production 



The inclusive production of the photon with a large transverse momentum in the 7p 
or ep collision at HERA and other colliders was the subject of detailed studies in the 
literature [8], [12]- [26], and herein only the main aspects of the calculation of the cross 
section including NLO QCD corrections are briefly discussed. The main cross section 
formula and explicit expressions for various contributions are given in Sec. 3.1 and in 
the Appendices. 

3.1 Cross section formulae 

The differential cross section for the ■jp —>■ jX (or jp —>■ ■y jet X) process has the 
following form: 

^^7P-7(iet)x ^ ^ /• dxfg^^ix, + / dxfg/p{x, iJ?)da^^-^^^ (3.1) 

q,q •' •' 

+ E [dx, Y f dxU/,{x,,iJi^)h/,{x,iJi^)da'^'-^' (3.2) 

a=g,g,9 b=q,q,g 

+ E fdx Y f^fwpi^,l^')D,/c{z,l^')da^'-'^'' (3.3) 

b=q,q,g c=q,g,g ■' ^ 

+ E / C^^7 E j dx Yl j ^fahi^'y,IJ''^)fb/p{x,IJ>^) 



dz 

0'=q,q,9 b=q,q,g c=q,q,g 



D^Jz,fi^)da'''^''' (3.4) 



+ j:fdx {fq/Ax, ^i') [cia^r"' + daZ-^^<^] + f,/,{x, ^i')daZ-^^^} , (3.5) 
g,g '' 

where 7g — > 75, — > etc. are various partonic processes described in the previous 
chapter. The functions /a/7, /b/p, and D^jf. stand for the parton distributions in the 
photon and proton, and the parton-to-photon fragmentation function, respectively. 
The corresponding longitudinal-momentum fractions carried by partons in the proton 
and photon are denoted as x and x-y. Similarly, for the parton-to-photon fragmentation 
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the variable z is introduced. The \i scale is the factorization/renormalization scale 
related to the physical hard scale of the considered process. 

The contribution of processes with the direct both initial and final photon are given 
on the right-hand side in (3.1) (see Figs. 2.1, 2.11). The expressions (3.2-3.4) stand 
for contributions of processes with the resolved photons or/and the fragmentation into 
the photon (Figs. 2.7-2.10). The formulae for the partonic cross sections, da^^^'^'' etc., 
appearing in (3.1-3.4) are listed in Appendix B. The 0{ols) corrections to the Born 
process are included in (3.5), where fic^l'^'''*, da']^^~^'^^ and da'^^'^'^^ are the partonic 
cross sections for the processes shown in Figs. 2.2-2.4. 

The 0{as) corrections to the Born process were first numerically calculated for the 
inclusive photon production in [15], and the compact analytical formulae can be found 
in [18, 25]. These corrections contain the infrared and/or coUinear singularities. The 
coUinear singularities are present in the squared amplitudes of the 2 — > 3 processes 
(Figs. 2.3, 2.4), and are factored out and absorbed by the corresponding parton den- 
sities in the proton or photon (Fig. 2.5), and by the parton-to-photon fragmentation 
functions (Fig. 2.6). The infrared singularities are due to the soft gluon emission or the 
soft gluon exchange and cancel when the real gluon corrections (Fig. 2.3) are added to 
the virtual gluon corrections (Fig. 2.2) in agreement with the Bloch-Nordsieck theorem 
[91] (see [18, 25] for the calculation details). The corresponding singular-free partonic 
cross sections da'^^^'^^ + da'2^^^^'^3 and da'^^^'^^'^ for the inclusive photon production 
are given in Eqs. (D.21) and (D.22) in Appendix D.2. 

To obtain the predictions for the inclusive photon production one needs to perform 
the integrations of the cross section (3.1-3.5) within the whole range of the fractional 
momenta: 0<x<l, 0<x^<l, 0<2;<1. The ©-functions in (D.21) and (D.22), 
and the 5-function in (B.l) ensure that in fact the x, and z close to zero do not 
contribute to the cross section. 

3.2 Numerical results 

The results for the inclusive (non-isolated) photon production at the HERA collider are 
obtained with the parameters and parton distributions given in Sec. 2.2. The initial 
electron and proton energies corresponding to energies at HERA are taken = 27.6 
GeV and Ep = 920 GeV, respectively [49]. 

In Fig. 3.1 the NLO cross section da/dE^, for the inclusive photon production in 
the ep scattering is presented (dashed line). The contribution of the processes with 
the resolved initial photon or/and the parton-to-photon fragmentation (solid line) and 
the contribution of the processes with the direct both photons (dotted line) are shown 
as well. The cross section decreases by three orders of magnitude when E'j. increases 
from 4 GeV to 20 GeV, and obviously the most important contribution to the total 
cross section is coming from the lowest Ej, region. The processes with the direct initial 
and final photons, i.e. Born plus 0{as) corrections and the box contribution (Figs. 
2.1-2.4, 2.11), dominate in the cross section for large Ej^ > 9 GeV. For lower Ej^ the 
processes with the resolved initial photon or/and the fragmentation into final photon 
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Figure 3.1: The cross section da/dEj, for the inclusive 7 production at the HERA coUider 
(dashed Une). The contributions of processes with the direct initial and final 7 (dotted line) 
and processes with the resolved initial 7 or /and the fragmentation into final 7 (solid line) 
are shown separately. 



(Figs. 2.7-2.10) constitute more than a half of the cross section. 

The importance of individual contributions to the inclusive photon cross section, 
integrated over 5 GeV < < 10 GeV, is illustrated in Tab. 1 (the first row) and 
in Figs. 3.2 and 3.3, where the distributions of the final photon rapidity, r]'^, are 
shown. In this E]^ range the integrated ("total") NLO cross section is equal to 241 pb. 
Processes other than the lowest order (Born) one give all together the contribution 
almost two times larger than the cross section for the Born process alone. The 0{as) 
corrections to the Born process are relatively small, negative at rj^ < and positive at 

> (Fig. 3.2), and constitute only 2% of the total cross section (Tab. 1). The box 
contribution constitutes 6% of the cross section integrated over all rj^ (Tab. 1) and it 
is even larger, 7-12%, for the differential cross section in the range — 1 < 77-y < 1 (Fig. 
3.2), where the data [42, 49] exists (see Chapters 4, 5, 6). Although the box diagram, 
7(7 — 7(7, is of order 0{a'^a1), it gives relatively large contribution mainly due to a 
large gluonic content of the proton at small x. 

The contribution of processes of order 0{a^a^) with the resolved initial photon and 
the quark-to-photon fragmentation is large, being of order 20% (Tab. 1). It is worth 
noticing that this contribution is build from many, relatively small, individual terms 
(Figs. 2.9, 2.10). The parton-to-photon fragmentation with the direct or resolved 
initial photon (Figs. 2.7-2.10) gives in sum a contribution of 31% (Tab. 1) or 32- 
39% for the central rapidity range — 1 < 77^ < 1 (Fig. 3.3), and this contribution 
is similar to the Born one (36%). Note that this part of the cross section involving 
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Figure 3.2: The cross section da/drf for the inclusive 7 production at the HERA coUider 
with (dashed hne) and without (dashed-dotted hne) the box contribution. The contributions 
of the Born process (sohd hne) and the Born process with the 0{as) corrections (dotted 
hne) are shown. 
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Figure 3.3: The cross section da/drp for the inclusive 7 production at the HERA collider 

(dashed line) and the contributions of various types of processes with: the direct initial 7 
and the fragmentation into the hnal 7 (short-dashed line), the resolved initial 7 and the 
fragmentation into the hnal 7 (solid line), the resolved initial and the direct hnal 7 (dashed- 
dotted line), and direct both 7's (dotted line). 
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poorly known fragmentation functions, D^^ is much suppressed by implementing of 
experimental cuts, see Chapter 4. 

In the whole range of r]^ the processes with the resolved initial and direct final 
photons constitute about 25% of the cross section ^ (Tab. 1, Fig. 3.3). It makes the 
cross section sensitive to the photon structure function [15], [18]- [24], [26]- [36], [92, 93] 
and in particular to the gluonic content of the photon, what will be shortly discussed 
in Chapter 7. 



^This contribution is even enhanced if experimental constrains are taken into account (Chapter 

4). 
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Chapter 4 

Isolated photon production. Small 
cone approximation 

In order to compare QCD predictions with data, one should consider such physical 
quantities, which are actually measured or are as close to the measured ones as pos- 
sible. Since in experimental analyses in order to reduce backgrounds the observed 
photon is often required to be isolated from hadrons [40]- [49], similar isolation of the 
photon produced in the ep scattering was included in the QCD calculations. It was 
done for the first time by Gordon and Vogelsang [27, 31] and Gordon [31]. The sec- 
ond independent calculation was performed by Krawczyk and Zembrzuski [28, 30]. 
Then the calculation of Fontannaz, Guillet, Heinrich, [29, 32] and a new analysis of 
Zembrzuski and Krawczyk [33] were presented. In this chapter we discuss the QCD 
predictions for the isolated photon production presented in the paper of Krawczyk and 
Zembrzuski [30] ^. 

4.1 Isolation restrictions 

The final photon is "isolated" if the sum of hadronic transverse energy within a cone 
of radius R around the photon is smaller than the photon transverse energy multiplied 
by a small parameter e [40] -[49]: 

El^"^"^ < eBTj.. (4.1) 

hadrons 

The cone is defined in the rapidity and azimuthal angle phase space, namely 

^ ^^hadron _ ^7)2 + {^(ffhadron _ ^7)2 < ^ (4.2) 

The isolation constraint allows to reduce the background from neutral mesons (tt^, rj), 
which decay into two photons and from photons radiated by final state hadrons. It 
suppresses considerably the contribution of processes involving the parton fragmenta- 
tion into the photon. For comparison with data the same isolation restriction (4.1) is 
taken into account in theoretical calculations [2 7] -[36]. 

^This is an extended and updated version of [28]. 
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Note, that a different type of isolation wliicli removes from tlie cross section tlie 
whole fragmentation component is advocated by Prixione [94]. We do not apply this 
type of isolation, since so far it is not applied in experimental analyses of the prompt 
photon production at the HERA coUider. 



4.2 Small cone approximation 

The simplest way to calculate the differential cross section for an isolated photon 
production, daisoi, is to calculate the difference of the differential cross section for a 
non-isolated (inclusive) photon production, danon-isoi, and a subtraction term, dasub 
[95, 96, 97, 27]: 

daisoi = danon-isoi — dasub- (4.3) 

The calculation of danon-isoi was discussed in the previous chapter. The subtraction 
term corresponds to the cuts which arc opposite to the isolation cuts, i.e. within a 
cone of radius R around the final photon the total transverse energy of hadrons should 
be higher than the photon transverse energy multiplied by e, 

^ ^hadron ^ ^^7_ (4 4) 

hadrons 

In this chapter we apply the subtraction term dasub calculated in an approximate 
way. This approximation bases on the assumption that an angle 5 between the final 
photon and a parton inside the cone of radius R is small [96, 27]. It allows to sim- 
plify the calculations considerably and leads to the compact analytical expressions for 
all relevant matrix elements involved in dasub- Note that in this approximation the 
maximal value of the angle 6 is proportional to the radius R: 

5 < R/cosh{r]^) = RE^/E-^, (4.5) 

where rj^, and Ej, stand for the photon rapidity, energy and transverse energy, 
respectively. We stress that the above small angle (cone) approximation is used only 
in calculation of the 0{as) corrections in the subtraction cross section dasub, while 
other contributions to dasub ^ well as danon-isoi are obtained in an exact way. 



4.3 Analytical results 

The calculations include partonic processes with two (2 — > 2) or three (2 3) particles 
in the final state: 

ab cd (4.6) 



or 



ab — 'ydid2, (4.7) 
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where a is the photon or a parton originating from the photon, 6 is a parton from the 
proton, c stand for the final photon or a parton which decays into the photon in the 
fragmentation process, and denotes quarks and/or gluons. For the isolated final 
photon the summation in Eqs. (4.1, 4.4) runs over the c-parton remnant and over the 

dj-partons, if they are inside the cone (4.2). 

In the partonic 2 — > 2 processes with a direct final photon (0 = 7), 

ab jd, (4.8) 

the photon is isolated by definition, so they give no contribution to the subtraction 
term. If the final photon comes from the fragmentation process (c 7^ 7), 

ab cd, (4.9) 

it takes the 2;-fraction of the c-parton transverse energy, Ej, = zE^. The hadronic 
remnant of the c-parton takes the fraction of transverse energy equal to (1 — z)Efj., 
and the photon is isolated if 

E:j^ = zE^ > e{l - z)E^ (4.10) 



or 



z>l/{l + e). (4.11) 

The subtraction term for 2 — > 2 processes can be obtained from cross sections (3.3, 3.4) 
involving fragmentation functions integrated over z, with z- values fulfilling inequality 
opposite to (4.11), i.e. z < 1/(1 + e). So we have 

b=q,q,9 c=q,q,g q 

1 1 ^'^^^^'^rf^ 

+ J2 12 dx J2 / —fah{x^,ll^)fblp{x,n'^) 

0'=q,q,g Q h=q,q,g Q c=q,q,g ^ ^ 

■D^/c{z,tJ,'')da'''^^. (4.12) 
The total subtraction cross section consists of da'j^'^g j^^ and 0{as) corrections: 

doZV'"" - daJZlL + . (4.13) 

The virtual gluon corrections do not contribute to the subtraction, since they come 
from 2 — > 2 processes (Fig. 2.2), where the direct final photon is isolated by definition. 
The cross section da^J'j^^^ includes only contributions of 2 ^ 3 processes (4.7). In 
these processes the photon and two partons arc produced. One parton, say di, enters 
the cone of radius R around the photon (4.2), and its transverse energy should be 
higher than the photon's transverse energy multiplied by e (4.4). There are three 
types of such processes, namely: 

•7? -^11 + 9 (with a quark inside the cone), (4-14) 
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•7? ~^ 19 + Q (with a gluon inside the cone), (4-15) 

•jg ■yq + q (with a quark or antiquark inside the cone). (4-16) 
The subtraction cross section for 0{as) corrections has a general form: 



V + vw 



U/,ix, /i^rfair^^-''^ + /,/,(x, /.^)rfaLV"''>-«^J dx, (4.17) 

where the variables v and w are defined in Appendix A. 2. These variables are related to 
the final photon momentum and to the initial momenta, and depend on the fractional 
momenta y and x. To fulfill the condition (4.4) the integration is performed over x 
according to © function with 

v{l-w)/{l-v + vw)> e, (4.18) 

since in the small cone approximation the transverse energy of the parton di is related 
to the photon transverse energy in the following way: 

E^' ^ E:I,v{1-w)/{1-v + vw). (4.19) 

The analytical results for the 0{as) corrections contributing to the subtraction 
term (4.17) are presented below. All coUinear singularities are shifted into the frag- 
mentation functions -D^/c (Fig. 2.6) according to the standard factorization procedure, 
as discussed in Introduction (for details for the non-isolated photon production see 
also [18, 25]). The infrared singularities do not appear in this calculations. This is 
due to the fact that there are no virtual gluon corrections, the gluons emitted inside 
the cone can not be too soft being constrained by Eqs. (4.18, 4.19), while the gluon 
emitted outside the isolation cone can not be soft, because its transverse momentum 
has to balance transverse momenta of the photon and the parton moving parallel to 
the photon. The cross sections for partonic processes (4.14-4.16), denoted in (4.17) by 
^^suV'^'^'^'^^ J are given by following expressions: 

E^^fl^ = e(I - .)g^<C. '\-"^7';t'^-f . P, (4.20, 

d-^p^ TTS^ [1 — V + vwy[l — V) 

d^p-y ^ ^ ns^ s 

(1 — v)\(l — V + vwY + (vw)'^],^ .4 4.^ .4, /, ^..x 

• (1 - . + v^Ml - ^) (L)^ ^ [1 + (1 - + + - - (4-21) 

d^p-r -^^^ 7r52 2 (l-^; + t;^)t;^(l-t;) ^' ^^"^^^ 
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where Cp = 4/3 and 



1 + - w) 



1 — V -\- vw 




+ 1. 



(4.23) 



The function P (4.23) which appears in Eqs. (4.20) and (4.22) was previously presented 
e.g. in [96], where it was used in calculations for the prompt photon production in 
the hadron-hadron collisions. Other terms on the right-hand side in Eqs. (4.20) and 
(4.22) correspond to the partonic cross sections for the 2^2 processes, 7g qg or 
19 ~^ QQj respectively. The expression (4.21) has been derived by us and presented 
for the first time in the paper of Krawczyk and Zembrzuski [30]. The derivation of 
formulae (4.20-4.22) is briefly discussed in Appendices D.5 and D.6. 

These results are obtained with the assumption that the angle between the photon 
and the parton inside the cone is small and only the leading terms arc kept. In 
d^a^r''"'^'^ (4.20) and ciV2y^'^'+^"' (4.22) the leading contribution is 0{lnR) + const 
and all terms of order 0{R^) or higher are neglected (terms of order 0{R) do not 
appear). In cPa]^l^"'^~^''' (4.21) there are no logarithmic 0{lnR) or constant terms. For 
small e, a soft gluon can be emitted into the isolation cone in the process 19 + Qi, 

leading to a large 0{R^) term in (4.21), and therefore this contribution can not be 
neglected. 

The small cone approximation was previously used to obtain the predictions for 
the prompt photon production at the hadron-hadron [96] as well as at the electron- 
proton [27] reactions. Then the exact (i.e. without the assumption that the cone is 
small) calculations were performed for both the hadron-hadron [98] and the electron- 
proton [31] collisions. The author of [98] compared the approximated and exact results 
and concluded that the applied approximation was reliable for the hadronic processes. 
Moreover, he mentioned in [31] that it is also an accurate technique for including 
isolation effects in the cross section for the prompt photon production in the electron- 
proton scattering. We have performed a detailed comparison between the predictions 
obtained for this process with the approximated and with the exact implementation 
of the isolation restrictions, see Sec. 5.4.1. 

4.4 Numerical results and discussion 

We perform the numerical calculations for the isolated photon production at the ep 

scattering using the HERA collider energies: Ee=27.5 GeV and Ep=820 GeV [42], 
and applying the parton densities specified in Sec. 2.2. In this section we use the 
parameter Aqc'£)=0.365, 0.320 and 0.220 GeV, fitted by us to the experimental value 
of as{Mz) = 0.1177 [99], for A^/=3, 4 and 5, respectively. 

4.4.1 Effects of isolation 

The results for the isolated photon production are presented in Figs. 4.1 and 4.2 in 
comparison with the results for the non-isolated photon production in the ep collision. 
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Figure 4.1: The differential cross section do/drp for the prompt 7 photoproduction at 
HERA as a function of the photon rapidity r]'^ for non-isolated photon (dashed line) and 
isolated photon with R = 1 and e = 0.05 (dotted line), 0.1 (solid line) and 0.2 (dashed- 
dotted line). The photon transverse energy, Ej, = pr, is taken in the range 5 < Ej, < 10 
GeV. The factorization/renormalization scale ji = Q = Ej, is used. 
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Figure 4.2: As in Fig. 4.1 for non-isolated photon (dashed line) and isolated photon with 
e = 0.1 and R = 0.1 (dashed-dotted line), 0.5 (dotted line) and 1 (solid line). 



The presented differential cross sections da/drf are obtained for 5 < E'j. < 10 GeV 
using various isolation cone parameters, e and R. The isolation cut suppresses the cross 
section by above 10% in the whole range of the photon rapidities shown in Figs. 4.1, 
4.2. For e=0.1 and i? = 1 the suppression is 17-23% at rapidities -1.5 < 77"^ < 4. This 
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large effect is not very sensitive to the value of e: changing the value by a factor of 2 
from e = 0.1 to e = 0.2 or to e = 0.05 varies the results for isolated photon by about 
4% (Fig. 4.1). The dependence on R is stronger but also not very large: when the 
R value is changed by a factor of 2 (from 1 to 0.5) the results increase by about 7% 
(Fig. 4.2). 

The suppression due to the isolation imposed on the photon is seen in Tab. 2, 
for the individual contributions and for the sum of all contribution. The results for 
both the non-isolated and isolated photons are shown in the first and the second row, 
respectively. As expected, the cross section for processes with fragmentation into final 
photon is strongly suppressed: due to the isolation it is lowered by a factor of 5. At 
the same time the 0{as) corrections to the Born diagram increase from 4.8 pb for the 
non-isolated 7 to 13.1pb for the isolated 7, i.e. the contribution of this corrections to 
the subtraction cross section, dasub, is negative. The requirement of isolation does not 
modify contributions of other processes since they are of 2 — > 2 type and involve only 
the direct final photons. The subtraction cross section, being a sum of negative 0{as) 
corrections and fragmentation contributions, is of course positive. Finally, the total 
cross section for isolated final photon is lower, by 20%, than for the non-isolated one. 

In further calculations discussed in this work the values R — 1 and e = 0.1 are 
assumed, as in analyses of the HI and ZEUS Collaborations [40]-[49]. 

4.4.2 Effects of other cuts 

In order to compare the results obtained by us with the data we take into account the 
isolation restrictions, as well as other kinematic cuts imposed in experimental analyses 
of the prompt photon events at the HERA collider. In this section we study effects 
of cuts applied in the analysis of the ZEUS Collaboration [42], where the initial 7 
fractional momentum and the final 7 transverse energy and rapidity arc required to 
be in the range: 0.2 < y < 0.9, 5 < Ej< < 10 GeV and —0.7 < rf < 0.9, respectively. 
The influence of the hmited y range is shown in Fig. 4.3. The cross section (sohd line) 
is strongly reduced, by 30-85%, in the positive rapidity region. At negative rapidities 
the change due to the y-cwi is weaker: 5-10% at —1.2 < f]"' < —0.4 and 10-30% at 
other negative rapidities. Separately we show the results obtained without including 
the box subprocess (dashed line). The box diagram contributes mainly in the rapidity 
region between -1 and 3. After imposing the y-cut it is important in the narrower 
region, from -1 to 1. The influence of the y-cut can be read also from the third row 
in Tab. 2: the Born contribution is reduced 3.5 times in comparison with the results 
obtained without the y-cut (the second row), while other contributions are suppressed 
less, roughly by a factor of 2. 

The results obtained for the isolated photon with the y-cut and in addition with 
the cut on the flnal photon rapidity, —0.7 < rj'^ < 0.9, are presented in the last 
row of Tab. 2. The restriction on t]'^ decreases the contributions of all subprocesses 
approximately by a factor of 2, except for the contribution of processes with resolved 
initial 7 and fragmentation into final 7 which is reduced almost 3 times. 

The influence of various cuts on the a;^-distribution is illustrated in Fig. 4.4. This 
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Figure 4.3: The differential cross section da/drf' for the photoproduction of the isolated 7 
(e = 0.1, R = 1) at HERA as a function of the photon rapidity rf with (solid lines) and 
without (dashed lines) the box contribution. The results are obtained with imposed y-cut 
(0.2 < y < 0.9) and without this cut. The photon transverse energy, Ej, = px, is assumed 
in the range 5 < < 10 GeV. 
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Figure 4.4: The cross section in bins of the length 0.1 for the photon transverse energy 
in the range 5 < Elj, < 10 GeV (Ej, = pt). The results for the non-isolated 7 integrated 
over the whole range of y and rf are shown with the dashed line. The solid line represents 
results integrated over the whole range of y and rf for the isolated 7 with e = 0.1 and 
R = 1. Results with additional cuts in the isolated 7 cross section are shown with: dotted 
line (0.2 <y<0.9) and dashed-dotted line (0.2 <y< 0.9, -0.7 < 77^ < 0.9). 
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distribution is scheme-dependent, nevertheless it is important to study the resolved 
photon contributions in the chosen renormalization scheme. Such an analysis was 
performed for the first time in our paper [30] . Note that Xj < 1 and Xj — 1 correspond 
to processes with the resolved or direct initial photon, respectively. So, the cross 
section in the last bin in Fig. 4.4, 0.9 < < 1, contain the contributions of both 
types of processes, while the cross sections in other bins, x^ < 0.9, consists of the 
resolved photon contribution alone. In particular we see in Fig. 4.4 that the isolation 
and the energy cut reduce considerably the cross section at large and medium Xj, 
while the cross section at x^ below 0.1 is reduced less. On the other hand, the small 
x^ contributions are strongly, by two orders of magnitude, diminished by the photon 
rapidity cut. This shows that measurements in the central rj'^ region (—0.7 < rj'^ < 0.9) 
are not too sensitive to the parton densities in the photon at small x^ region. 

4.4.3 Dependence on the choice of the renormahzation scale 

In order to estimate the contribution due to missing higher order terms, the influence of 
the choice of the rcnormalization/factorization scale /i is studied for the r]'^ distribution. 
In Fig. 4.5 the results obtained with and without the y-cut are shown. When changing 
H from E'j. to 2-E'j. {Ej^/2) the cross section increases (decreases) at rapidities below 
~ 1 and decreases (increases) at higher rapidity values. Only at high rapidities, rj'^ > 3, 
the dependence on a choice of the scale is strong, above 10%, however here the cross 
section is relatively small. In the wide kinematic range, —2 < rj'^ < 2, the relative 
differences between the results obtained (with and without the y-cut) for n = Ej, and 
lor jjL — 2 ■ Ej, OT Ej,/2 are small and do not exceed 6%. Around the maximum of 
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Figure 4.5: The cross section do/drp for 5 < Ej, = pT < 10 GeV and fi = Q = Ej,/2 
(dashed hne), Ej, (sohd hne) and 2 ■ Ej, (dotted hne). The results obtained without and 
with the y-cut, 0.2 <y < 0.9, are shown. 
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the cross section at rapidities — 1 < r;'^ < these differences are 4-6%. This smaii 
sensitivity of the results to the change of the scale is important since it may indicate 
that the contribution from neglected NNLO and higher order terms is not significant. 
This is important for a comparison (see below, and in Sec. 5.5.3 and 6.2.4) with other 
NLO calculations using other set of diagrams and, in particular, including the 0{as) 
corrections to the processes with the resolved photon which are not included in our 
calculation. 

Note that individual contributions are strongly dependent on the choice of /x, e.g. 
results for the processes with the resolved photon vary by ±10-20% at rapidities rj'^ < 1 
while the 0{as) corrections to the Born diagram vary by a factor of a few and in some 
kinematic ranges they even change sign (not shown). The results are much more stable 
when the sum of the resolved photon processes and the 0{as) corrections to the Born 
process is included (compare the discussion in Sec. 2.1.2). 



4.4.4 Comparison with the ZEUS data 

The results of the first measurements of the isolated prompt photon cross section 
(with no jet requirements) in photoproduction events at the HERA collider have been 
published in [42] . Below we compare of our predictions with these data. 



Fig. 4.6a shows the distribution of the photon transverse energy, Ej,, for various num- 
ber of active massless quarks, Nf. The predictions for Nf = 4 and 5 are in agreement 
with most of the experimental points but tend to lie slightly below the data. The 
difference between predictions for Nf = 4 and Nf = 5 is very small due to the fourth 
power of electric charge characterizing processes with two photons. The contribution 
of the charm quark is an order of magnitude larger than the contribution of beauty and 
the predictions for Nf = 3 are significantly below the predictions for Nf — 4 and below 
the data, confirming the need to include the charm contribution in the calculations. 

77^, Nf 

A similar comparison between the predictions and the data, now for the distribution 
of the photon rapidity, rj^, is shown in Fig. 4.6b. A good description of the data is 
obtained for Nf=A, 5 in the rapidity region 0.1 < rf < 0.9. For —0.7 < rf < 0.1 
our predictions lie mostly below the experimental points. Note that this disagreement 
between predicted and measured cross sections is also observed for calculations of 
Gordon [31] and of Fontannaz, Guillet and Heinrich [29] as well as for Monte Carlo 
simulations [42, 45]. In order to find the source of such a disagreement various checks 
have been performed, discussed below. 
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Figure 4.6: Tiie cross section da/dEj, (a) and da/drp (b) with various numbers of active 
massless flavors: Nf = 3 (dashed hnes), 4 (soUd hnes) and 5 (dotted hnes), compared to 



the ZEUS data [42]. 
(dot-dashed hne) 



The result without the box contribution is also shown for Nf = 4 



rf^, integrated a, C(a^a|) 

First, we study the size of 0{a^a'g) contributions included by us. In Fig. 4.6b we 
present separately an effect due to the box subprocess (for Nf — 4). One can see that 
the box term enhances considerably the cross section in the measured rapidity region. 
Its contribution to the integrated cross section is equal to 10%, see Tab. 2 (fourth 
line). The processes with resolved photon and fragmentation also give a no- negligible 
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contribution, although roughly two times smaller than the box one (Tab. 2, fourth 
line). We conclude, that taking into account of these 0{a'^al) contributions improves 
the description of the data. 

Next, we study the dependence on the parton densities. The predictions obtained 
using three different NLO parton densities in the photon (ACFGP^ [24], GRV^ [72] 
and GS^' [81]) are presented for A^/ = 4 in Fig. 4.7a = E'^) and in Fig. 4.7b 
(/i = 2E:1-) together with the ZEUS data [42]. The results based on ACFGP'^ and 
GRV''' paramctrizations differ by less than 4% at rapidities 7]"' < 1 (at higher t]'^ the 
differences are bigger), and both give a good description of the data in the rapidity 
range 0.1 < 77^ < 0.9 (for = and = 2EJ). For -0.7 < 77^ < 0.1 none of the 
predictions is in agreement with the measured cross section. 

For fi=Ej' (Fig. 4.7a) the GS''' distribution leads to results considerably below ones 
obtained using ACFGP^ and GRV^ densities, especially in the rapidity region from 
roughly -1 to 1, where the differences are up to 14%. This difference between the GS''' 
and other considered herein parton paramctrizations is mainly due to their different 
treatment of the charm quark in the photon, namely in the GS''' approach the charm 
quark is absent for /i^ below 50 GcV^. Since we take 5 < /i — Ej^ < 10 GeV, and 
the most important contribution to the cross section arises from the lower Ej^ region 
(see Fig. 4.6a), the fj,"^ value usually lies below the GS''' charm quark threshold. As 
a consequence, the predictions based on GS'"' have strongly suppressed contribution 
of subprocesses involving charm from the photon - contrary to GRV' and ACFGP'"' 
predictions where the charm threshold is assumed at lower /x^. 

All the considered parton paramctrizations give similar description of the data 
when the scale is changed to /x = 2Ej., see Fig. 4.7b. Here the calculation corresponds 
to /x^ which is always above 50 GeV^ and the charm density in the GS''' parametrization 
is non-zero, as in other paramctrizations. 

V'^, fp, 

We have also studied the dependence of our results on a choice of the parton densities 
in the proton and parton fragmentation into the photon (not shown). The MRST1998^ 
(set ftOSa) [74] and CTEQ4M^' [77] NLO densities in the proton give predictions lower 
than GRV^ [71] NLO densities by 4-7% and 4-6% at negative rapidities. At positive 
rapidities the cross sections calculated using MRST1998^, CTEQ4]VP and GRV^ vary 
among one another by less than 4%. Results for the isolated photon production are 
not too sensitive to the parton-to-photon fragmentation function. For rapidity ranging 
from -1 to 4 the cross section obtained with DO-^''"^(LO) [15] fragmentation function 
is 2 — 4% lower than the cross section based on GRV'^^'^^(NLO) [73] parametrization. 
Only at minimal [r]'^ < —1) and maximal (4 < rj'^) rapidity values this difference 
is larger, being at a level of 4 — 8%. All these parton densities in the proton and 
fragmentation functions lead to similar descriptions of the data. 
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Figure 4.7: The cross section da/drp in the range 5 < Ej, < 10 GeV for the renormal- 

ization/ factorization scale fi = Ej, (a) and fi = 2 ■ E^ (b). The GRV^ [72] (solid hne), 
ACFGP^ [24] (dotted line) and GS^ [81] (dashed line) parton densities in the photon were 
used together with GRW densities in the proton [71] and GRV^'''"'^ fragmentation [73]. 



rj'^, y 

In Fig. 4.8 our predictions are compared to the ZEUS data divided into three ranges 
of the initial photon fractional momentum y. This allows to establish that the above 
discussed discrepancy between the data and the predictions for r]'^ < 0.1 is coming 
mainly from the low, 0.2 < y < 0.32, and medium, 3.2 < y < 0.5, y region. In the 
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Figure 4.8: The cross section da/dri"' for three ranges of y: 0.2 < y < 0.32 (a), 0.32 < y < 0.5 
(h) and 0.5 < y < 0.9 (c). The predictions of K&Z [30] and LG [31] are compared with the 
ZEUS Collaboration data [42]. The figure is taken from [42]. 

high y region, 0.5 < y < 0.9, a good agreement is obtained. 



summary of comparison with data 

To summarize, there are discrepancies between the predictions and the ZEUS data [42] 
at r]'^ < 0.1 and y < 0.5. The variation of theoretical parameters and parton densities 
do not improve the description of the data. Note that there are also discrepancies 
between predictions and the new HI data for the prompt photon production at HERA 
[49] (Chapter 5). So far the source of these discrepancies is not established. 

Finally, it is worth mention that in Ref. [43] there was performed an analysis, in 
which our computer program for the prompt photon production at HERA [30] was 
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modified: the parton densities in the photon for a;^ > 0.8 were multiplied by a factor 
of 2 or 3. For doubled at > 0.8 the description of the data was much improved, 
and for tripled /y at > 0.8 a good agreement with the ZEUS data was found in the 
whole range of rj"'. This agreement is probably accidental, but it can also indicate a 
need of improvement of the parton densities in the photon at Xj close to 1 [43] . 

4.4.5 Comparison with other QCD predictions (LG) 

The calculation of L. Gordon (LG) [31] for the photoproduction of isolated photons 
at HERA contains the 0{as) corrections to the processes with the resolved initial 
photon, which are not included in our calculation (Chapter 2). On the other hand, 
we include the box diagram (Fig. 2.11) neglected in [31]. The 0(a|) diagrams with 
resolved initial 7 and fragmentation into final 7, as well as other diagrams shown in 
Figs. 2.1-2.4 and 2.7-2.10, are taken into account in both calculation. 

The LG calculation apphes the MRST^(NLO) [74] parton distributions in the pro- 
ton which give predictions 4% lower on average than GRV^(NLO) [71] parametrization 
used by us (Sec. 4.4.4). Both calculations use GRV'*' [72] and GS'*' [81] parton densities 
in the photon. Since the GS''' densities lead to worse description of the data than the 
GRV' ones (see Sec. 4.4.4), below we discuss the predictions obtained using the GRV' 
parametrization. 

The comparison between the predictions of Krawczyk and Zembrzuski (K&Z) [30] , 
predictions of Gordon (LG) [31] and the data was performed by the ZEUS Collabora- 
tion [42]. Results of this analysis are shown in Figs. 4.8, 4.9 and 4.10 ^. 

The LG predictions for da/dEj, cross section are higher than ours in the whole pre- 
sented range of transverse momentum, 5 < Ej, < 15 GeV (Fig. 4.9). Both calculations 
agree with the data within large experimental errors, although the K&Z predictions 
tend to he shghtly below the data. 

For da/drj^ the biggest differences arc at ?7'*'=0.9 where the LG results are about 
35% higher than K&Z (Fig. 4.10b). The differences decrease towards negative rapidity 
values and are negligible at —0.7 < rf < —0.5. The average difference in the range 
—0.7 < rf < 0.9 is about 10%. Both predictions are in disagreement with data at 
rf < 0.1. At rj'^ > 0.1 K&Z agrees with the measured cross section but LG tends to 
lie shghtly to high. 

For y range limited to low values only, 0.2 < y < 0.32, the LG cross section is 
higher than ours by up to 20% at positive rj'^, while at negative r]"^ it is lower by up to 
10%. (Fig. 4.8a). In this y range, as well as for medium range with 0.32 < y < 0.5 
(Fig. 4.8b), the K&Z and LG calculations lead to similar discrepancies with data at 
rj'^ ^ 0.1. For large y values, 0.5 < y < 0.9, where our predictions agree with data, the 
LG results are higher than ours by up to 70% at 0.7 < rf < 0.9 and does not describe 
the data well (Fig. 4.8c). 

^In the comparison between LG [31] and K&Z [30] we consider the predictions presented in the 
ZEUS paper [42], because the numerical results presented in the original paper [31] are given for 
different kinematic ranges. 
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Figure 4.9: The cross section da/dEj, for the photoproduction of isolated photons in the 
kinematic range 0.2 < y < 0.9 and -0.7 < rf< < 0.9. The predictions of K&Z [30] and 
LG [31] using GRV^ [72] parton distributions in the photon, are compared with the ZEUS 
Collaboration data [42]. The Monte Carlo simulations are also shown. The figure is taken 
from [42]. 
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Figure 4.10: As in Fig. 4.9 for the da/drj^ cross section in the range 5 < < 10 GeV. 
In addition, the predictions of K&Z [30] and LG [31] using GS^ [81] parton densities in the 
photon are shown. The figure is taken from [42]. 



39 



CHAPTER 4. ISOLATED PHOTON PRODUCTION. SMALL CONE 
APPROXIMATION 



To summarize, the differences between both calculations including different sets 
of diagrams are of order 10%. However there are some kinematic ranges where the 
differences are larger up to 70% for 0.5 < y < 0.9 and 0.7 < rf < 0.9. The exper- 
imental errors arc to large to conclude which calculation describe better the ZEUS 
Collaboration data presented in [42]. Nevertheless, for large y region, 0.5 < y < 0.9, 
our predictions give better description of the data. 

There is also another NLO calculation for the photoproduction of isolated photons 
at HERA presented by Fontannaz, Guillet and Heinrich (FGH) [29] . We compare our 
predictions with predictions of FGH in Sees. 5.4.3 and 6.2.4, where we use the same 
parton distributions as FGH do. Let us only note here, that in the range considered 
by the ZEUS group [42], 0.2 < y < 0.9 and -0.7 < t]^' < 0.9, the FGH predictions 
obtained using MRST1999^' [75], AFG^ [79] and BFG^'''^^ [83] parametrizations are 
about 10% higher on average than our predictions obtained using GRV^''"''-^''"^ [71, 72, 
73] parametrizations. However, for the large y region, 0.5 < y < 0.9, and for r)'^ in the 
narrowed range 0.7 < 7]"' < 0.9 the FGH predictions are about 50% higher. 

The dependence of our predictions on the renormalization/factorization scale n is 
about ±5% in the kinematic range considered by the ZEUS Collaboration [42] and 
it is true also for 0.5 < y < 0.9. So, we could expect that the size of higher order 
contributions is not large in this particular region of y as well. However we found that 
the difference between K&Z and LG (FGH) is up to 70% (50%) for 0.5 < ?/ < 0.9 and 
0.7 < T]"' < 0.9. To understand the origin of this large difference we should know the 
dependence of the LG and FGH predictions on the scale /x in this particular range, 
which however was not presented in [31] and [29] ^. If this dependence was large, it 
could be a confirmation of our supposition that the 0{as) corrections to the processes 
with resolved photons should be included together with the 0(a|) corrections to the 
Born process in ncxt-to-ncxt-to-lcading order (NNLO) accuracy and, that they should 
not be included in NLO to avoid possible uncertainties due to the choice of ii. 



^The dependence on the choice of the scale ji was studied in [31, 29] for the total cross sections, 
i.e. cross sections integrated within considered therein kinematic ranges. The total cross section of 
FGH for the prompt photon production vary by about ±4%. The total cross section of LG for the 
prompt photon plus jet production vary by ±3%. 
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Isolated photon production. Exact 
implementation of isolation cuts 

In the small cone approximation, discussed in the previous section, one obtains simple 
analytical formulae which are easy to implement in numerical calculations for the 
isolated photon production. The accuracy of this approximation is good [96, 27, 98, 31], 
however the method has one drawback: it does not allow to obtain any predictions 
for the rapidity or transverse energy of the jet (or jets) produced together with the 
photon. This is because the formulae (4.20-4.23) and (D.21,D.22) used in the small 
cone approximation are integrated over all allowed momenta of final particles other 
than the photon. 

In order to compare predictions with the existing data not only for the photon 
observed in the final state but also for the ep — > 67 jet X process, there is a need of 
more general technique allowing to consider both the photon and the jet momenta in 
any kinematic range for which the data are taken. An introduction of such a technique, 
based on the division {slicing) of the three-body phase space [100]-[107], [29, 31, 32] 
is presented below and was discussed in the paper of Zembrzuski and Krawczyk [33]. 
Note, that our method of the division of the phase space [33] shghtly differs from 
methods applied in other papers. 

The use of the phase space slicing allows to calculate the cross section in an exact 
way, i.e. with no assumption that the isolation cone is small. We take some cross 
section and matrix elements formulae from the literature and derive all the missing 
expressions, see Sees. 5.1-5.3 and Appendices B-D. The comparison between exact 
predictions for the ep — > e^X cross section and the predictions based on the small cone 
approximation, as well as the comparison with the new HI data [49] for the ep — > e'jX 
reaction, are presented in Sec. 5.4. 

The slicing of the phase space is applied to obtain predictions for the ep ej jet X 
cross section, and comparing them with corresponding HI data [49], in Chapter 6. 
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5.1 Cross section formulae 

The general formula for the isolated photon (or isolated photon + jet) cross section is 
obtained from Eqs. (3.1-3.5) integrated over the whole range of x and and over z 
within the region given by Eq. (4.11), namely: 

1 1 



d^iP^i{3et)x ^Y^l dxf^/^{x, t/)da'^'i^^'^ + / dxfg,^{x, fi^)da^^^'^^ (5.1) 

1 1 

+ E jdx.Y. J dxfa/^{x^,fi')fb/p{x,fi')da'''^'''' (5.2) 



a-q,q,g q b=q,q,g g 

1 1 



+ Y fdx Y: I %h/Ax,lj')D,/,{z,y?)da'^'--'^ (5.3) 
b=g,g,g o ^=9'9'%/i+e ^ 



dz 



+ Y [ dxj Y f dx Y I -^fah{x^^l^'^)h/p{x,iJ? 

a=q,q,gQ b=q,q,g q c=q,q,g^^^^^ ^ 



■D^,,{z,n')da'''^"' (5.4) 

1 

+ Yldx {/,/p(x, //^) [d<7^™ + d^™] + /,/,(x, /.^)d<--^^^"} , (5.5) 

q,q 

where the partonic cross sections in (5.1-5.4) are given by Eqs. (B.1-B.22). Note, that 
these partonic cross sections are not integrated over the photon momentum, so they 
allow to obtain various types of differential cross sections (e.g. da/drp or da/dEj^) in 
any range of the final photon or final parton momenta. This is due to the fact that in 
2 — > 2 processes the parton momentum is determined, in the centre-of-mass system, 
by the final photon momentum. The isolation is taken into account simply through 
the lower limit of the integration over z. 

An inclusion of the isolation in calculations of the 0{as) corrections (5.5) is more 
complicated, since one needs to restrict the final momenta of two partons, and to take 
care of a cancellation and factorization of singularities. 



5.2 Division of phase space 

The 0{as) corrections (5.5) include 2 — 3 processes (4.7) and the virtual gluon cor- 
rections. The squared matrix elements of these processes contain infrared and coUinear 
singularities [18, 25], see Appendix D. In order to obtain numerical predictions for the 
isolated photon production with arbitrary kinematic cuts, one needs to isolate these 
singularities, and cancel or factorize them in such a way, that an integration over any 
range of final momenta is still possible [29, 31]. The standard method to achieve this 
goal is to divide the phase space (Appendix C) into a few parts [100]-[107], [29, 31, 32], 
and herein this technique is adopted. 
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Let us assume that 6ji is an angle between the momentum of the final parton di 
(see Eq. (4.7)) and the momentum of an initial particle j, where j = e, p ^. Let us 
also define the distance, R^i, between the di parton and the final photon: 



We define the parts of the phase space in the following way: 

• Part 1. In the first part the variable w (see Appendix A.l) is in the range 
uJcut < < 1, where Wcut is an arbitrary parameter close to 1: < 1 — Wcut ^ 1- 
This region of the phase space contains all types of NLO corrections: the virtual 
gluon exchange, the real gluon emission, and the process ^yg — > ^yqq. The virtual 
gluon exchange is a 2 — > 2 process, while the other ones are of 2 — 3. However, 
for w close to 1, Wcut < w < 1, the two final partons in the 2^3 processes are 
almost coUinear or/and one of the final partons is soft. For Wcut sufficiently close 
to 1 the kinematics of the 2 — > 3 processes is almost the same as in the 2 — > 2 
case: 



where the "particle" d has four-momentum equal to the sum of four-momenta of 
the dj-partons, and its mass is almost zero. The additional restriction 1—Wcut < e 
is sufficient to ensure that if parton di or d2 is inside the isolation cone then its 
energy is smaller than eE^, so the final photon is isolated. 

The other parts of the phase space. Parts 2-5 described below, are defined for w < Wcut, 
and contain only 2 — > 3 diagrams. 

• Part 2. In this region w < Wcut and 6cut > niin{9ei^9e2)i where Ocut is a small 
arbitrary parameter, Ocut ^ 1- Here, one of the final partons has the momentum 
almost coUinear to the momentum of the initial electron and it, for sufficiently 
small Ocut) does not enter the isolation cone around the photon. The second 
parton has a large transverse momentum balancing the photon transverse mo- 
mentum, so the final photon is isolated. 

• Part 3. Here w < Wcut ^-nd Ocut > min{0pi,0p2). One of the final partons has 
momentum almost coUinear to the momentum of the proton, and, like in Part 
2, the final photon is isolated. 

• Part 4. Here w < Wcut a-nd Rcut > min{Rji, R-y2), where Rcut is a small param- 
eter, Rcut -C 1. In this case one of the final partons, say di, is almost coUinear 
to the photon, and the photon is isolated if the transverse momenta obey the 
inequality (4.1): Exdi < eEx-y 

^The initial particles a and b in Eq. (4.7) move parallel to the initial electron and proton, 
respectively. 




(5.6) 



76 70? 



d = di + d2, 



(5.7) 
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• Part 5. The last part is defined as the region in which there are no colhnear 
configurations: w < Wcut, Ocut < min{9ei,9e2), O^ut < min{9pi,9p2), and Rcut < 

The 0{as) corrections (5.5) contain contributions from all parts of the phase space: 

'daZr^' + daZ-'^'] = [dall-f + daH^r] + E ^^^^ (^-S) 

1=2,3,4,5 

rf^7|->7..-^ E ^<^r (5-9) 

1=1,2,3,4,5 

In numerical calculations for Part 1 the formulae (D.21) and (D.22) (the same as for the 
inclusive photon case) can be used with the requirement w > Wcut- In these formulae 
all the soft gluon singularities present in the virtual gluon and real gluon corrections 
are canceled, and all the colhnear singularities are factorized into the parton densities. 

Parts 2-5 contain no contribution from the virtual gluon exchange, and no soft 
gluon singularities. Analytical expressions for cross sections corresponding to Parts 2, 
3, and 4 are derived in Appendix D, and will be presented in next section. 

The cross section in Part 5 has no singularities at all, and one can perform an 
exact numerical integration over final four-momenta in any kinematic range, including 
isolation restrictions and other cuts. For this purpose the general formulae (C.13) and 
(D.1-D.20) with £ = are used. 

The presented method of the division of the phase space differs from the methods 
applied in other calculations for the photoproduction of isolated photons at HERA 
[31, 29, 32], since we use the cut-off parameter Wcut not used before. Our method 
allows to obtain relatively the simplest analytical singular-free formulae in all parts of 
the phase space, as all the singularities due to the soft gluons (or gluons colhnear to 
the final quark) are collected in one part of the phase space, i.e. in Part 1, and this 
part is chosen this way that the formulae for the non-isolated photon production can 
be used. The unintegrated over momenta formulae used in Part 5 are also the same 
as for the non-isolated photon case, and in Parts 2-4 the relatively simple expressions 
for colhnear configurations are applied. 

The drawback of our method is that the numerical calculations are not stable if 
the integration over final momenta in Part 5 is performed in the laboratory ep frame 
or in the rest frame of the initial photon and parton. However, we found that the 
calculations were sufficiently stable if the rest frame of the final partons was chosen 
instead. 

The predictions should not depend on a choice of unphysical cut-off parameters, 
1 —Wcut, 9 cut and Rcut, if they are small enough. On the other hand, they can not be too 
small, since for very low values of these parameters large numerical errors occur. The 
results presented in Sees. 5.4 and 6.2 are obtained with 1 — Wcut = (^cut = Rcut = 0.01 
where the angle 9cut is defined in the centre of mass of the initial photon and the initial 
parton (despite the fact that the numerical calculation is performed in other frame). 
We have checked that the change of the predictions due to the variation of these 
parameters is negligible, below 1%, if they are taken in the range 10^'^ < 1 — Wcut ^ 
0.03, 10"^ < 9cut < 0.05 and 10"^ < Rcut < 0.05. 
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5.3 Analytical results 

Below we present analytical results for cross sections for 2 — > 3 processes including 
various collinear configurations described in Sec. 5.2. These results are obtained (for 
the details see Appendices C, D) using the standard factorization procedure, which for 
the inclusive (non-isolated) photon production in the ep scattering was discussed in 
Ref. [18, 25]. All the collinear singularities appearing in the calculations are factorized 
into parton densities in the photon (Part 2), proton (Part 3), and into fragmentation 
functions (Part 4). 

Some of the presented below "splitting functions" P have been already presented in 
Ref. [96] , where they were used to calculate in the small cone approximation the cross 
section for the isolated photon production in the hadron-hadron collision (see also e.g. 
[101, 103, 104]). Nevertheless, we find it useful to present herein compact analytical 
formulae for all the relevant cross sections for the isolated photon production in the 
ep scattering. The partonic cross sections for the 2^3 processes in the parts of the 
phase space labeled in Sec. 5.2 as Part 2, 3 and 4 have the form: 

1 

d^'sS'' = e{s + t + u)l da'i^-^^\^s, ^t, u)P,^,g{^, yEe{l - O0cut)d^, (5.10) 



1 

dc^asS"' = e{s + t + u)J da'^^-'^^^s, ^t, u)P,^,q{i, yEe{l - i)e,ut)di, (5.11) 



1 

^C^r -^{s + t + u)j da-^i^^^^s, t, e«)P™(e, xEp{l - O0cut)d^, (5.12) 



1 

d^^Zir = e{s + t + u)J da^'^-'T.^^s, t, Cu)Pg^,g{i, xE,{l - O0cut)d^, (5.13) 



t + u 



.e(s + t + u) |d<j™(s,-,|)P,^^,(e,£;?(l-eRn*)^, (5.15) 



t + u 



Q{s + t + u) J da^^-^\s,^-,^)P,^,,(^,EUl - ORcut)f, (5.17) 



re 
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where + e) includes the isolation restrictions (4.1), the partonic cross sections 

for the 2 — > 2 processes are given in Appendix B and the functions P are given by 



W«.^J = ^4{ii^lnf (5.18) 



P^^,,(e, E^) = 3— e le + (1 - 01 In ^ + 1 , (5.19) 



P,^,,(e, E^) = ^^^^^(1 - E^), (5.20) 



^.-^..-(e, - (5-21) 



The variable ^ in the functions Pa-^bc{^: E_i) stands for the fraction of the energy of the 
particle a taken by the particle b: E^ = ^Ea- The variable E±_ is the maximal possible 
(for given dcut or Rcut) momentum of the particle c perpendicular to the momentum 
of the particle a: E± = maxlpc x Pa|/|Pa|- The results (5.10-5.21) are obtained with 
the assumption that 9cut ^-nd Rcut are small, see Appendices C and D for details. 

Note, that the integration over ^ in Eqs. (5.10-5.17) is trivial as the 2 — > 2 partonic 
cross sections, da"-^^^'^{s,t,u), contain the S{s + t + u) function, sec Eq. (B.l). The 
expressions (5.15) and (5.17) are equivalent to the expressions (4.20) and (4.22) for 
the subtraction terms used in the small cone approximation (Sec. 4.2). 

The cut-off parameter Rcut is chosen in a Lor cntz- invariant form (Sec. 5.2), but 
the parameter 9cut depends on the choice of the frame of reference. It can be defined in 
any frame of reference provided that the energies yEe and xEp in Eqs. (5.10-5.13) are 
given in the same frame. In numerical calculations we use the expressions (5.10-5.13) 
with 6cut defined in the centre of mass system of the initial photon and the initial 
parton, where yEg = xEp. 



5.4 Numerical results and discussion 

The predictions for the photoproduction of the prompt photons at the HERA collider 
are obtained with parameters and parton densities specified in Sec. 2.2. For a com- 
parison with the HI Collaboration data, kinematic ranges used in [49] are applied: 
the fraction of the electron energy transferred to the photon is restricted to the range 
0.2 < y < 0.7, the final photon rapidity and transverse energy are taken in the limits 
— 1 < Tj'^ < 0.9 and 5 < Ej. < 10 GeV, and the initial energies are Eg = 27.6 GeV and 
Ep = 920 GeV. The final photon is isolated with parameters e = 0.1 and R—1. 
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5.4.1 Comparison between exact and approximated results 

In Figs. 5.1 and 5.2 the differential cross sections da/dEj. and dcr/drj'^ are presented. 
The exact predictions ^ (sohd hnes) are compared with the predictions obtained in 
the small cone approximation (dashed lines) and with the new HI data [49] (see also 
[48]). The differences between both predictions are about 3%, so the small cone 
approximation works well, despite the fact that the isolation cone of radius R — 1 is 
not a small one. 

The comparison between the two predictions are also presented in Tab. 3 for 
individual contributions and total cross sections. The Born, box and resolved photon 
contributions are the same in both calculations, since the approximation was applied 
to 0{as) corrections only (Sec. 4.2). Nevertheless the difference between exact and 
approximated 0{as) corrections is very large. We have found that for the isolated 
photon production the 0{as) corrections to the Born process obtained using the small 
cone approximation are 67% higher than the corresponding contribution obtained 
using the division of the phase space ^. We conclude that the small cone approximation 
gives a very rough estimation of the order of magnitude of the 0{as) corrections, 
however the total cross section is estimated quite well, since the 0{as) corrections 
give relatively small (5%) contribution to the total cross section (Tab. 3). 

5.4.2 Comparison with the HI data 

The comparison with the new HI Collaboration data [49] is presented in Figs. 5.1 
and 5.2 and we see that the predictions tend to lie below the data for r]"' distribution, 
although in most bins an agreement is obtained within the large experimental errors 
(Fig. 5.2). For Ej. distribution (Fig. 5.1) the data in two bins are 2-3 standard 
deviations above the predictions, while in other bins the the data and predictions 
agree well. 

Our predictions shown in Figs. 5.1, 5.2 are obtained for the final photon and 
partons production (Eqs. 4.6, 4.7). According to the analysis done by the HI Collab- 
oration [49], for a realistic comparison with their data, the hadronization corrections 
and multiple interactions (h.c.-|-m.i.) should be taken into account in the predictions. 
Our exact predictions, after h.c.+m.i. corrections implemented by the HI Collabora- 
tion [49], are shown in Fig. 5.3 (dotted lines) together with the HI data [49] and with 
predictions of Fontannaz, Guillet and Heinrich (FGH) [29, 32]. Both, our and FGH, 
predictions are obtained with with MRSTQQf [75], AFC [79], BFG-'^'^"^ [83] parton 
densities. To show the effect of h.c.-|-m.i., the FGH predictions are presented with 
(solid lines) and without (dashed lines) these corrections 

The FGH and our predictions are below the data typically by 30% and 40%, re- 
spectively [49]. The discrepancy is seen in almost whole range of 1]"' with exception 

^By exact predictions we mean the predictions obtained using the technique described in Sec. 5.2 
with no assumption that the isolation cone is small. 

^For non- isolated photon there is no difFcrcnce between both types of calculations. 

^The correction factors are typically 0.75-0.90 depending on a bin. In both K&Z and FGH pre- 
dictions the same factors are applied. 
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Figure 5.1: The comparison between exact predictions (solid line), predictions obtained in 
the small cone approximation (dashed line) and the HI [49] data for da/dEj, cross section. 
The GRVP^'^'f''"-^ [71, 72, 73] parton densities are used 
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Figure 5.2: As in Fig. 5.1 for da/drj^ cross section. 
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Figure 5.3: The cross sections da/dEj- (a) and da/drj"' (b) for the ep — *■ e'jX process 
with -1 < r/T < 0.9, 5 < < 10 GeV and 0.2 < |/ < 0.7. The HI data [49] are com- 
pared with K&Z [33] (dotted hnes) and FGH [29, 32] (sohd hnes) predictions obtained 
using MRST99P [75], AFG^ [79], BFG-''™^ [83] parton densities. Both predictions are 
corrected for hadronization and muhiple interactions (h.c.+m.i.) effects. The shaded 
areas show changes of FGH predictions due to the fi scale variation between 0.5-£'j' 
2 ■ Ej^ (inner bands) and added hnearly uncertainties on h.c.+m.i. corrections (outer 
bands). The LO (dashed-dotted) and NLO (dashed) FGH results without h.c.+m.i. 
are also shown. The figure is taken from the HI paper [49] 



of the last bin, where the QCD calculations and data are in agreement. Note, that 
the prompt photon cross sections measured previously by the ZEUS Collaboration 
[42] were well described by predictions at 77''' > 0.1 (Sec. 4.4.4). On the other hand, 
both the ZEUS [42] and the HI [49] data are consistent, as it was shown in [49, 48]. 
It seems, that the above different conclusions of comparisons between the predictions 
and the ZEUS data, and between the predictions and the HI data are due to large 
effects of h.c.+m.i. corrections applied in the HI analysis [49]. 

Note, that the dependence on the choice of parton parametrizations is relatively 
large: results obtained using GRV^''^'-^'"'^^ [71, 72, 73] (not shown) are 10% higher on 
average than these obtained with MRST99P [75], AFG^ [79], BFG-^'^'^f [83] densities, 
and slightly decrease the size of discrepancies between the predictions and the data. 
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5.4.3 Comparison with other QCD predictions (FGH) 

The calculation of Fontannaz, Guillet and Heinrich (FGH) [29, 32] for the photopro- 
duction of isolated photons at the HERA collider includes the 0{as) corrections to 
the processes with the resolved initial photon or/and with the fragmentation into the 
final photon, which are not included in our calculation (Chapter 2). Both our and 
FGH calculations take into account the (9(a^a|) box diagram (Fig. 2.11) and all 
other diagrams shown in Figs. 2.1-2.4 and 2.7-2.11. 

The predictions of FGH presented in their papers [29, 32] were obtained in the 
kinematic ranges considered by the ZEUS Collaboration [42, 44] (see Sec. 4.4.5) and 
one can not use them for the direct comparison with our predictions corresponding 
to the HI measurements [49]. So, to compare the predictions, we use the FGH and 
K&Z (Zembrzuski and Krawczyk) [33] results based on the theoretical predictions and 
corrected for h.c+m.i., as presented in the paper of the HI Collaboration [49] ^. Both 
predictions are obtained using the same parton parametrizations, namely MRST99^ 
[75], AFGT [79], BFG-^'-"^ [83]. 

The FGH predictions are 15-20% larger than ours for the r)"' and Ej, distributions 
shown in Fig. 5.3. and this difference is of a similar order as the difference between 
our predictions and the predictions of Gordon (LG) [31] which was discussed in Sec. 
4.4.5. It means, that the higher order terms included in the LG and FGH calculations 
and not included in the K&Z approach, give a sizable contribution to the cross section 
in the kinematic range considered for the prompt photon photoproduction. This is not 
fully supported by our study of the dependence on the renormalization/factorization 
scale, which seemed to suggest smaller effect due to missing higher order terms (Sec. 
4.4.3). 



^The applied correction factors are the same for both K&Z and FGH predictions, see footnote ^. So 
the relative differences between both "corrected" predictions are the same as the differences between 
pure QCD (with no h.c.-|-m.i.) K&Z and FGH predictions 
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Chapter 6 

Isolated photon plus jet production 



The partons produced in hard processes (Figs. 2.1-2 A, 2.7-2.11) are not observed 
experimentally, since they recombine into colorless jets. In previous parts of the work, 
the predictions for the prompt photon production were discussed and the presented 
cross sections were integrated over partons/jets momenta. Now, the isolated photon 
plus jet production is studied for limited ranges of both photon and jet rapidities and 
transverse energies. Such a process with two objects observed in the final state can be 
a source of more detailed information concerning the dynamics of the interaction. 

6.1 Jet algorithm 

As it was discussed in Chapter 2 and in Sec. 4.2, in the calculation we deal with 2 — > 2 
and 2 — > 3 hard processes. In 2 — > 2 processes (4.6) the jet originates from the parton 
d and the momentum of the jet can be identified with the momentum of this parton. 

Two partons, di and ^2, produced in 2 — 3 processes (4.7) may lead to two separate 
jets or they may form one jet. The number of observed jets and their momenta 
depend on the jet definition. In this work the inclusive /cy-jet finding algorithm [108, 
109] is appfied in accordance with the jet definition used in the HI Collaboration 
measurements [49]. 

Since in the presented calculation there are no more than two partons forming the 
jet or jets, the algorithm becomes very simple. If the distance between the partons, 
i?i2, defined as 



is larger than an arbitrary parameter then two separate jets arise with transverse 
energies, rapidities, and azimuthal angles of the dj-partons: 




(6.1) 




n 



Jet; 




I = 



1, 2. 



(6.2) 



For Ri2 < Rj both partons are treated as components of one jet with 



Jl/rp H/rp -\- H/rp , 



(6.3) 
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r^ie* = (E^^T)'^^ + E^^ri'^^ ) lEi^\ (6.4) 

= (£;^i(/,<ii + E^^(t)'^^)/E{^\ (6.5) 

The algorithm can be easily applied in numerical calculations using the phase space 
slicing (Sec. 5.2). In the first part of the phase space, where w ~ 1, the kinematics is 
the same as in 2 —> 2 processes. In part 2 (3, 4) one parton moves parallel to the initial 
electron (initial proton, final photon) and it, for sufficiently small Ocut {Ocut, Rcut), does 
not enter the cone defining the jet, so the jet consists (on the partonic level) of the 
second parton alone. In part 5 one or two jets may arise depending on the value of 
Ri2 (6.1). 

Following experimental analyses [49], i?j = 1 is used in our numerical calculations. 



6.2 Numerical results and discussion 

We perform numerical calculations in kinematic regions as used in experimental anal- 
ysis. There are two publications of the ZEUS Collaboration presenting results of 
measurements of the isolated photon plus jet photoproduction at the HERA col- 
lider [40, 44]. In the first paper [40] the cross section integrated over some kinematic 
range is given. The aim of the second one [44] was to study transverse momentum of 
partons in the proton, and no results for cross sections were presented (distributions 
of events, not corrected for the detector effects, were shown). There is also a confer- 
ence paper of the ZEUS Collaboration where data for a differential cross section for 
the isolated photon plus a jet photoproduction are given, however these data are still 
preliminary [41]. 

In the new paper of the HI Collaboration [49] (see also [48]), the photoproduction 
data for the isolated final photon (e = 0.1 and R = 1) with the initial energies 
Ee — 27.6 GeV and Ep = 920 GeV are presented for various differential cross sections 
of both ep — > ejX (considered in Chapter 5) and ep — > 67 jet X processes. Herein, 
the predictions for the kinematic limits as in [49] are compared with these recent final 
data. The cross sections are integrated over 0.2 < y < 0.7, and —1 < rf < 0.9 
and/or 5 < Ej^ < 10 GeV with the jet rapidity and jet transverse energy in the range 
— 1 < T]^^^ < 2.3 and 4.5 GeV < respectively. If two jets are found within the 
above region, that with higher E^* is taken. Other parameters are specified in Sec. 
2.2. 

6.2.1 Asymmetric cuts 

As it is discussed in [32], the symmetric cuts for the photon and the jet transverse 
energy, E^l^^^ — E'j.^j^, lead to unphysical results in next-to- leading or higher orders 
of calculations due to constraints imposed on soft gluons (see also [110, 111, 112]). 
We have decided to study this effect for the cross section as a function the photon 
transverse energy, not analyzed in [32]. This allows to understand the effect of the 
symmetric cuts in more details. The results of this study was first presented in our 
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paper [33]. Next, similar study of the effect of the symmetric cuts was presented 
in [34]. 

Our findings are illustrated in Fig. 6.1a,b. We study the dependence of the NLO 

predictions on the photon transverse energy in the i?^-range wider than the range 
considered by the HI Collaboration. At Ej, values close to the minimal jet transverse 
energy, = 4.5 GeV, the NLO differential cross section has a discontinuity, see 

Fig. 6.1a. For {Erp)_ — > 4.5 GeV the cross section has a strong maximum while a 
minimum for (Ej,)^ — > 4.5 GeV. In the minimum the value of the cross section is even 
negative. This unphysical fluctuation is due to processes with soft gluons: large terms 
corresponding to the soft gluon emission and the virtual gluon exchange do not cancel 
properly if the photon transverse energy is close to the minimal transverse energy of 
the jet, because some of this terms corresponds to Ej. slightly below the cut-off E^l^^^ 
and the other correspond to Ej^ slightly above Ej^l^^^. 

However, if we integrate the cross section over some E'j.-bins then this fluctuation 
disappears. As it is shown in Fig. 6.1b, the NLO predictions are well defined if one 
takes the cross section integrated over E^ from Eipl^-^^ — A to Eij^l^-^^ + A, provided 
that A is sufficiently large in the comparison with the gap in Fig. 6.1a, say A > 0.3 
GeV. So, the bins of a length 1 GeV presented in Fig. 6.1b are large enough to avoid 
errors corresponding to the fluctuation around Ej, — E^^^^. 

An integration of the differential cross section (Fig. 6.1a) over the photon transverse 
energy higher than the minimal jet transverse energy, E^ > E^l^^^ (symmetric cuts), 
leads to underestimated predictions in NLO. However numerically this effect is not 
very important, being at level of roughly 5% if the integration is performed in the 
range from E^ = = 4.5 GeV to E]. = 10 GeV. 

Note that in the previous measurements [40, 41, 44] the symmetric cuts Ej, > 5 
GeV and > 5 GeV were used, while in the new HI analysis [49] the asymmetric 
cuts E:}. > 5 GcV and E^^^^ > 4.5 arc taken. This latter choice of cuts allows to avoid 
the theoretical errors in QCD calculations and is applied in calculations presented in 
next sections (Sec. 6.2.2-6.2.4). 

It is worth mentioning that the cross section for Ej^ < E^^^^^ (Fig. 6.1a) is dom- 
inated by the contribution due to the 2 — 3 processes in the 0{as) corrections 
(Figs. 2.3, 2.4), since the contributions of the 2—^2 processes (Figs. 2.1, 2.7- 
2.11) are suppressed by the requirement E^^ > E^l^-^ and by the isolation. For 
Ej, < £'^'^^j^/(l + e) 4 GeV the 2^2 processes do not contribute at aU. On the 
other hand, for larger E'jf*, say Ej^^ > 5 GeV, the contribution of 0(as) corrections 
is very small, see Tab. 1. 



6.2.2 Theoretical uncertainties 

The dependence of the considered cross sections on the choice of the renormalization/ 
factorization scale, /i, is not strong: variations of fj, from Ej. to Ej,/2 or 2Ej, lead to 
changes of the cross section less than 3% for Ej, < E^l^^^ and up to 5% for Ej, > E^l^^^ 
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Figure 6.1: The differential cross section da/dEj, (a) and the differential cross section 
da /dEj, averaged over Ej,-hins (b) for the ep — > ej jet X process. The NLO predictions for 
Ej, = n (a) and for jj, between Ej,/2 and 2Ej, (b) are shown. 



(Fig. 6.1b). The predictions for various are also shown in Fig. 6.2. In each presented 
here bin the dependence on the choice of fi is less than ±5% for both (Fig. 6.2a) 
and Tjjet (Fig. 6.2b) distributions. The total cross section integrated over the kinematic 
range considered by the HI Collaboration [49] vary by ± 3.4%. Since the effect of the 
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Figure 6.2: The cross section da/dE^^ (a) and da/drf^^ (b) for /x = 0.5 • Ej. (dashed-dotted 
lines), n = Ej, (solid lines) and /x = 2 • Ej, (dashed lines). The HI Collaboration data [49] 
are shown for comparison. 



variation of the ii scale is not large, the calculation seems to be stable, and one can 
expect that the contributions of higher orders are not sizable. 
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Figure 6.3: The cross section da/d'q'^ for MRST2002P (NLO) [76] (dotted line), GRW 
(NLO) [71] (solid line) and CTEQ6P (NLO) [78] (dashed line) parton densities in the proton 
used with GRV^'f^'^^ (NLO) [72, 73] densities in the photon and fragmentation functions. 
The HI Collaboration data [49] are also shown. 



u 

Next, we have checked the sensitivity to the choice of parametrizations. In Fig. 6.3 

the results obtained using different parton densities in the proton arc shown. The 
predictions of CTEQ6P (NLO) [78] are 6% higher than the predictions of MRST2002P 
(NLO) [76]. The GRV^ (NLO) [71] densities give results higher than MRST2002P by 
5-7% at negative r/^, and 3-5% at positive r]^. Diffcrencies between CTEQ6^ and GRV^ 
do not exceed 4%. 



The most important variables for testing the structure of colliding particles are the 
fractional momenta of partons in these particles. Below we consider the distribution 
of the fractional momentum in the photon, however since the theoretical variable 
is not a good observable, in experimental analyses some estimations of are used 
instead. We consider here the observable x'^^ , which is defined as [40, 47] ^: 

x°^' = {Ei;'e-'>'^' + E].e-''')/2yEe. (6.6) 

The variable equal to the "theoretical" one, x^, for 2 ^ 2 processes with the direct final 

photon. For the processes with a larger number of partons in the final state and for the processes 
with the parton-to-photon fragmentation, the x^ and x°''^ differ. 
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Figure 6.4: Tiic cross section da/dx^"' for GRW (NLO) [72] (solid line), AFGOT/ (NLO) [80j 
(dashed line) and OS"' [81] (dotted line) parton densities in the photon used with GRW'^'^"'^ 
(NLO) [71, 73] densities in the proton and fragmentation functions. 



In Fig. 6.4 the x°°^ distributions are shown for different parton densities in the photon. 
The GS^ (NLO) parametrization [81] give predictions lower than GRV^ (NLO) [72] 
by 20-36% for < 0.9. This large difference is due to the specific treatment of the 
charm contribution in the GS parametrization (see Sec. 4.4.4). The resuhs obtained 
with use of AFG^ (NLO) [79] and AFG02^ (NLO) [80] are very similar, and only 
the latter is shown in Fig. 6.4. It gives predictions up to 15% lower than GRV' for 
x°^^ < 0.9. At l&rge-x"^^ region, x°^^ > 0.9, the cross section is dominated by the 
contribution of processes with the direct initial photons, and the differences between 
predictions obtained using various parametrizations are small. For the total cross 
section integrated over all x^^^, within the considered range 
the difference between results based on GRV^ and AFG02^ (G^S) is 4% (16% ). 

Note, that there are new LO [82, 113, 114] and NLO [115] parametrizations of 
the real photon structure with a special treatment of heavy quark contributions. For 
the first time for the photon they use the ACOT^ scheme introduced previously for 
the proton [116] and include the newest F2 data, never used in constructing other 
parametrizations for the photon. We have compared the results obtained using the 
LO parametrization of Cornet, Jankowski, Krawczyk and Lorca (CJKL'*^) [82] with 
results obtained using GRV'''(LO) [72]. Despite the fact that parton densities in both 
parametrizations differ considerably, they lead to similar results for the prompt pho- 
ton production at HERA. Predictions for CJKL'''(LO) are about 3% lower than the 
predictions for GRV'''(LO) (in this comparison the 0{as) corrections were not taken 
into account and the GRV^(LO) [71] densities in the proton were used). 
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We have also compared predictions of DO^^'^^ (LO) [15], GRV-'^'^"^ (NLO) [73] and 
BFG-'^''"^ (NLO) [83] fragmentation functions. The isolation requirement reduces the 
contribution of processes involving the resolved final photon (Sec. 4.4.1, Tabs. 1, 
2), so the dependence on the choice of fragmentation functions is weak, even if the 
fragmentation functions differ considerably from one another. The total cross sections 
for the isolated photon plus jet production obtained with DO'^^"^ and BFG-''''"^ (set I 
and set II) are lower than the predictions of GRV-^''"^ by 2% and 4%, respectively. 

fp^ f-y^ 

The GRV distributions for the proton [71], photon [72] and fragmentation [73] have 
been used as a reference in calculations discussed above, and while performing the 
comparison each time only one paramctrization has been changed. The differences 
observed in the total cross section (i.e. in the cross section integrated within the 
considered kinematic range) are not large (with an exception of the GS densities, 
which give predictions considerably lower than the other densities in the photon, as we 
discussed before). However, the differences can be larger if one changes simultaneously 
all the used distributions. For instance predictions of the MRST1999^ [75], AFG''' [79] 
and BFG-^™f [83] set ^ arc lower than the GRVP-^'-^^^prcdictions by 10% on average. 
The comparison between both results is presented in Fig. 6.5. The differences are 
large, up to 13%, for GeV, while for higher Ej, both sets of parton densities 

give predictions relatively close to each other and differences are below 6%. 

In the calculation the renormalization/factorization scale fi = E"^ is used for Ej^ be- 
tween 5 and 10 GcV (with exception of Fig. 6.1, where a wider rage of Ej^ is shown). 
It is a standard assumption that the number of active massless flavours at this scale 
is Nf — A [27]- [34]. For a comparison, in Fig. 6.6 the results obtained with Nf — 3 
and Nf = 5 are also shown. The contribution of the bottom quark is very small 
and the results for Nf = 4 and Nf = 5 are similar. Differences between these two 
results are much smaller than the standard deviations of the HI data [49] which are 
also presented in Fig. 6.6. On the other hand, the contribution of the charm quark 
is large: the results obtained using Nf — 3 are about 35% below the predictions for 
Nf — A. Neglecting of the charm mass may lead to a slight overestimation of the 
production rate, especially in the box contribution which is particularly sensitive to 
the change from Nf = 3 to Nf = 4. However we do not expect that an improved 
treatment of the charm contribution would change our results qualitatively, since the 
energy /j, — Ej, > 5 GeV is several times larger than the charm mass. 

^ These parametrizations are used for comparison with the HI data and with the FGH predictions 
in Sees. 5.4.2, 5.4.3, 6.2.3 and 6.2.4. 
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Figure 6.7: As in Fig. 5.3 for the ep ej jet X process. The additional cuts for the jet 
are: Si^^ > 4.5 GeV and -I < 77^'^* < 2.3. The figure is taken from [49]. 



6.2.3 Comparison with the HI data 

The QCD results for the photon plus jet production shown in Figs. 6.2, 6.3, 6.5 and 
6.6 are in reasonable agreement with the data of the HI Collaboration [49] although 
some discrepancies are present, especially for r]jet < —0.3 (Fig. 6.2) and for the Ej, 
distribution (Figs. 6.5, 6.6). 

However, according to [49], for the realistic comparison with the data, the pure 
perturbative QCD calculations should be corrected for effects of hadronization and 
multiple interactions (h.c.+m.i.), as discussed in Sec. 5.4.2. Such a comparison was 
performed by the HI Collaboration [49], see Figs. 6.7 and 6.8, where the predictions of 
Zembrzuski and Krawczyk (K&Z) [33] (dotted lines) corrected for h.c.+m.i. are pre- 
sented together with the HI Collaboration data [49] and with predictions of Fontannaz, 
Quillet and Heinrich (FGH) [29, 32]. To show the size of the corrections, the FGH 
NLO results are plotted without (dashed lines) and with (solid lines) h.c.+m.i ^. Both 
K&Z and FGH predictions are obtained using the MRST1999P [75], AFG^ [79] and 
BPQ/''as [83] parton parametrizations. 

The K&Z and FGH calculations for the ep —>■ 67 jet X process give somewhat 
better description of the data than in the case of the ep — > e'yX process (Sec. 5.4). 
The agreement is seen in most bins in Figs. 6.7 and 6.8, especially when one takes 
into account on the one hand the uncertainties of h.c.+m.i. corrections and on the 
other hand the theoretical uncertainties due to the variation of the fj, scale (±5%) and 



^See footnote in Sec. 5.4.2. 
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Figure 6.8: As in Fig. 6.7 for da /dE^r^'' (a) da/dr]^^^ (b) da/dx!^^ (c) and da/dx^^ (d) cross 
sections. The figure is taken from [49]. 



to the choice of parton densities (~ 10%) (Sec. 6.2.2). Nevertheless, the K&Z and 
FGH predictions tend to underestimate the HI data. In some kinematic ranges the 
data are 1-2 standard deviations above the K&Z predictions e.g. for 6.7 < Ej. < 9.2 
GeV (Fig. 6.7a) and 0.3 < r]jet < 1-6 (Fig. 6.8b). 
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In Figs. 6.8c,d the cross sections da/dx^^ and da/dx^'^ are shown, where x^'^ and 
x^'^ are defined as follows [49] ^: 



X. 



(6.7) 

(6.8) 



The K&Z predictions are 1-1.5 standard deviations below the data for x^^ < 0.85 (Fig. 
6.8c). Note, that better agreement is obtained if GRV^'''''-'^''"^ [71, 72, 73] parametriza- 
tions are used (not shown in Fig. 6.8 which is taken from the HI paper), since they 
give predictions 18% higher than MRST1999P, AFG^ and BFG^™^^ at x^^ < 0.85 (at 
> 0.85 the difference is 4%). 

At x^^ > 0.012 {x^^ < 0.0034) the K&Z predictions are 2 standard deviations 
below (above) the data (Fig. 6.8d). An implementation of GRV^'''''-^''"^ densities does 
not improve the description of the data (not shown). 

The HI Collaboration has also presented the cross sections da/dp± [49], where p± is 
the component of the photon momentum perpendicular to the scattering plane (see [49] 
for the precise definition). For p± the 2 — > 2 processes give no contribution, so 
the cross section is sensitive to higher order processes only. For x;^*^ > 0.85 the cross 
section is dominated by 0{as) corrections to the processes with direct initial and final 
photons, which are included in NLO in our calculation, and the data are in reasonable 
agreement with K&Z predictions. On the other hand, for x:^*^ < 0.85 the cross section 
is dominated by 0{as) corrections to the processes with resolved photons. These 
contributions are not included in the K&Z calculation, so our predictions by definition 
can not describe such data. 

6.2.4 Comparison with other QCD predictions (FGH) 

The calculation of Fontannaz, Guillet and Heinrich (FGH) [29, 32] takes into account 
the 0{as) corrections to the resolved photon processes, which are not included in 
the calculation presented in this work (Chapter 2, Sec. 5.4.3). In the considered 
kinematic range for the photon plus jet production the total cross section of FGH is 
about 4% higher than our predictions, so the total contribution of 0{as) corrections 
to the resolved photon processes is relatively small. 

The K&Z and FGH results shown in Fig. 6.7b differ by 5% or less in the whole range 
of rj'^. The differences are larger for other differential cross sections (Figs. 6.7a, 6.8). At 

> 6.7 GeV the FGH predictions are about 13%. higher, while at 5 < < 5.8 GeV 
they arc 7% lower than predictions of K&Z (Fig. 6.7a) The largest differences are 
for 1 < 1]^'^^ < 1.6 and 1.6 < 1]^'^^ < 2.3, where the FGH predictions are above K&Z 
by 33% and 63%, respectively. Despite these divergences, both calculations lead to a 
similar description of the HI data shown in Figs. 6.7 and 6.8. 

^For 2^2 processes the experimental variables xi^"^ and Xp"^ are equal to corresponding theo- 
retical variables multiplied by z (defined in Sec. 3.1): xi^'~'=z ■ Xj, Xp'^=z ■ Xp. 
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Note that the FGH predictions are closer to the HI data than the K&Z predictions 
for da/dp± with xi^'^ < 0.85 [49], since we do not include processes which contribute 
in this region (Sec. 6.2.3). 
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Chapter 7 

Probing the gluon content of the 
photon 

The photoproduction processes in the electron-proton scattering are sensitive to the 
parton densities in the proton as well as in the photon and provide an opportunity to 
probe the photon structure at the HERA collider (for a review of the data see [3]). 
Also the photoproduction of photons with a large transverse momentum (the Deep 
Inelastic Compton process) was considered as a possible tool to test the quark and 
gluon densities in the photon and proton [20]- [24], [26]- [33], [92, 93]. In particular, it 
was found that the contribution involving the gluon density in the photon is especially 
large when the final photon is produced in the forward (proton) direction [20]- [24], 
[92, 93]. Recently new analysis devoted to the possibility of measuring the gluon 
density in the photon and proton in the ep 67 jet X process at HERA were 
presented in [34, 35, 36]. 

So far we have considered the production of the photon with the large transverse 
momentum (transverse energy) in the electron-proton scattering where the mediating 
photons have been quasi-real, ~ 0, and they spectrum have been given by the 
equivalent photon approximation [86]-[89] (Sec. 2.3). Now, we extend the study 
including the effect of small, but non-zero virtuahty, 7^ 0. 

The first attempt to describe the Deep Inelastic Compton process at HERA using 
the parton densities in the virtual photon can be found in [117, 118], where wc tested 
in LO the validity of the equivalent photon approximation. Next, we have examined 
the usefulness of this process to study at the HERA collider the structure of a virtual 
photon, and in particular its gluonic content [92, 93]. In this chapter we discuss shortly 
some of our LO results obtained in [92, 93]. 
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7.1 Calculation of the cross section 

The invariant differential cross section for the deep inelastic electron-proton scattering 
can be written in the following form: 

£;:^^^ = r(a^>-- + eaZ>--), (7.1) 



where Uj. [al ^~*^) is the cross section for the interaction between the proton and 
the virtual photon polarized transversely (longitudinally) and p'^ {E'^) stands for the 
final electron momentum (energy). Coefficients F and eF are functions of energies and 
momenta of the electron in initial and final states, sec [9]. They can be interpreted 
as the probability of emitting by the initial electron the virtual photon polarized 
transversely and longitudinally. If so, we can use (7.1) to obtain the differential cross 
section for the ep — > e^X process taking into account the virtuality of the exchanged 
photon: 

^ ^^ep^ejX ^^7'p^7X ^^^*p^^x 

where Py (E^) stands for the final photon momentum (energy). Since the cross section 
for the reaction ep — > e'-fX is dominated by the exchange of photons with small 
virtuality, one can neglect a contribution due to the longitudinal polarization [5]. 
Assuming that the exchanged photon has only the transverse polarization we obtain: 

^^^^ ^ i ^^-^^ 



^7 



Our aim is to study the sensitivity of the cross section to the gluon distribution in 
the photon. To achieve this goal we include the Born process (Fig. 2.1) and the 
processes with the resolved initial photon (Fig. 2.7), and we omit contributions of 
other processes: 

^^^Fp " 5 i dxf,,p{x, + (7.4) 



■^7 

q,q 

1 



r f d(T°'^~*'^'^ 
dxy Yl / dxfa/y'{xy,n^,Q'^)fb/p{x,iJ?)Ey— , (7.5) 

o.=q,q,9 b=q,q,g ^ ^ Pi 

where f^jp and /a/7* are the parton densities in the proton and in the virtual photon. 
In the partonic cross section da^**^^"*^ as well as in the partonic distribution /a/7* we 
take the virtuality of the photon emitted by the electron exactly as it follows from 
the kinematics of the process. 
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7.2 Numerical results 

In the calculation we use the GRS (LO) [84] and the GRV (LO) [71] parton distribu- 
tions in the virtual photon and in the proton, respectively. The number of flavours 
is assumed Nf = 3 - the maximum number of active massless quarks in the GRS 
parametrization. We take the QCD parameter Aqcd = 0.2 GeV and the hard scale 
equal to the transverse energy of the final photon, n — Ej,. Calculations are performed 
for the initial electron and proton energies in the HERA accelerator: Ee = 30 GeV^ 
and Ep = 820 GeV^ with the virtualities of the mediating photon, ranging from 
10~^ GeV^ to 2.5 GeV^. Note that due to a smooth behaviour of the GRS parametriza- 
tion in the hmit ^ we were able to perform the calculation also for below 

We study the differential cross section 

for the transverse energy of the final photon E^ = 5 GeV and the energy of the 
mediating photon E^* = yEg with fixed y — 0.5. 

As it was expected, the cross section decreases for increasing values of the virtuahty 
Q^. For example the predictions obtained for = 0.25 GeV^ (2.5 GeV^) are one order 
(three orders) of magnitude lower than the predictions for = 0.03 GeV^ [92]. 

The processes initiated by the gluonic content of the virtual photon dominate in 
the cross section in the forward (proton) direction, i.e. for large photon rapidities, 
rjj [92, 93]. This effect is presented in Fig. 7.1, where the contributions due to the 
resolved virtual photon divided by the Born contribution are shown as a functions of 
rjy. The results for the process gytQp — * 79 are obtained taking between 10"'' GeV^ 
and 1 GeV^. For comparison we show the results for q^*gp — > jq and Q'7*^p(^7*?p) — jg 
processes obtained using = 0.1 GeV^. 

The large dominance of the contribution due to the gluonic content of the virtual 
photon over the Born contribution is seen at large 77, in the whole considered range of 
the photon virtualities. The largest ratio, up to factor of 35, is obtained for the small 
virtualities. It is important since in the limit — > the value of the cross section is 
the largest. Thus the cross section integrated over Q"^ between and, say, 1 GeV^ ^ 
is also strongly dominated by the process g^Qp — > 75 in the forward direction. 

The process involving the gluonic content of the photon dominates at large rj^ not 
only over the Born contribution but also over the contribution of processes involving 
quarks from the photon, see Fig. 7.1 and Ref. [92]. This gives an opportunity of 
measuring the gluon content of the real and virtual photon at HERA. 



^This range of is considered for the photoproduction of prompt photons at HERA, see Sec. 
2.3. 



66 



7.2. NUMERICAL RESULTS 
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Figure 7.1: The results for the ep e'-fX reaction including non-zero virtualities, Q^, 
of the mediating photon obtained for y=0.5 and E]. = 5 GeV. The curves represent 
the cross section (7.6) for the g'y*qp 7g (solid lines), q'y*gp — >■ jq (long-dashed line) 
as well as q^qp —>■ jg and q^qp 19 (short-dashed line) processes divided by the 
corresponding Born contribution. The predictions for g^*qp — > 7g process are obtained 
with = 10~^, 0.01, 0.25, 0.5 and 1 GeV^, while for the other processes involving 
resolved photons — 0.1 GeV^ was used. 
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Chapter 8 
Summary 



We have presented a NLO QCD calculation for the photoproduction of the isolated 
photons and the isolated photons associated with jets at the DESY HERA ep collider. 

Our calculation includes set of diagrams different than other existing NLO QCD 
calculations for this process. This difference arises from different treatment of parton 
densities in the photon and parton-to-photon fragmentation functions. In our approach 
these parton densities and fragmentation functions are of order 0{a) while in most of 
other calculations they are treated as quantities being of order 0{a/as). Our counting 
originate from the fact that the "structure" of the photon and the fragmentation into 
the photon arises from the purely electromagnetic processes j ^ qq and q — > jq, 
respectively. It ensures a cancellation to a large extend of dependences on the choice 
of the renormalization/factorization scale /i in the cross section. 

We have presented a method to divide (slice) the three-body phase space which 
differs from methods applied in other calculations for the considered process. Our 
method allows to obtain relatively simple analytical singular-free formulae for the 
0{as) corrections in each part of the phase space. Dividing of the phase space is 
introduced in order to implement exactly isolation restrictions and other kinematic 
cuts. It was used to obtain predictions for the isolated photon and isolated photon 
plus jet production at HERA. 

The numerical predictions obtained using the phase space slicing were compared 
with the predictions obtained using the small cone approximation. In this approxima- 
tion the isolation cuts in the 0{as) corrections to the Born process are implemented 
in an approximated way while in other contributions the isolation is implemented ex- 
actly. We confirm the observation of other author, that the small cone approximation 
gives very accurate predictions. However we have found that it is due to the fact that 
the 0{as) corrections to the Born process are small. These 0{as) corrections alone 
are obtained within the small cone approximation with very low accuracy. 

The theoretical uncertainties of our predictions are well under control. The predic- 
tions vary within 10% for various used parton parametrizations. The dependence on 
the choice of the renormalization/factorization scale fi for the isolated photon (isolated 
photon plus jet) productions is below ±6% (±5%) in each considered kinematic region, 
and the average dependence is ±5% (±3%). This small sensitivity to the choice of /i 
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may indicate that the not included contributions of higher orders (NNLO or higher) 
are not sizable. 

Other NLO QCD calculations (including some diagrams which are of NNLO order 

from our point of view) give predictions 10-20% (4%) higher on the average than ours 
for the isolated photon (isolated photon plus jet) production. However there are some 
kinematic ranges where the differences are larger, up to 70% for the isolated photon 
production in the range 0.5 < y < 0.9 and 0.7 <rf < 0.9 (where our predictions give 
the best description of the data) . This large difference is not supported by our study 
of the sensitivity to the choice of ^, which seems to suggest smaller effect of higher 
order contributions. To understand these differences we should know the dependence 
of other predictions on the scale ii in this particular range, which however was not 
provided. 

We have studied in details the influence of the isolation restrictions as well as other 
experimental cuts. In particular, we have shown that the measurements in the central 
rapidity range considered in current experiments are not much sensitive to the parton 
densities in the photon in the low-a;^ region. 

Effects of symmetric cuts for the photon and jet transverse energies were studied. 
We have found that the NLO differential cross section da/dE], is a discontinuous 
function at Ej, = Ej^^^^^. This is due to constraints imposed on the phase space of 
soft gluons. We confirm the known result that to avoid theoretical uncertainties the 
asymmetric cuts are preferred. 

Our predictions and computer program were used by experimental groups for com- 
parison with data. The QCD predictions for the photon plus jet production tend to 
lie below the data, nevertheless they agree with the data in most of bins. On the other 
hand, none of existing predictions describes the data satisfactory well for the photon 
production with no cuts for the jet. It means that the cuts for the observed jet remove 
from the cross section contributions of such kinematic configurations for which the 
disagreement with the data and the differences between predictions are the largest. 

Finally, we have pointed out the sensitivity of the Deep Inelastic Compton cross 
section to the gluon density in the photon also for the virtual photon. We presented a 
formula which allows to include the parton densities in the virtual photon in the cross 
section for the ep — > ejX (or ep — > ej jet X) reaction. We have found that, as it 
was expected, the processes initiated by the gluon arising from the photon dominate 
in the cross section in the forward direction, i.e. for large rj^. 
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Tables 



The tables contain our predictions for the ep — > e'yX and ep 67 jet X photopro- 
duction processes at the HERA colhder. The initial electron and proton energies, the 
final photon isolation parameters (e = 0.1, i? = 1) as well as the kinematic ranges of 
y, r/T, ri^^\ and E^f* are taken from the HI [49] (Tabs. 1, 3) and ZEUS [42] (Tab. 
2) Collaborations papers and are described in the captions of corresponding tables. 

The following notation is used for various contributing processes: 

• 0{as) = 0{as) corrections to the Born process, 

• dir,dir = direct initial and final 7, 

• dir,frag = direct initial 7 and fragmentation into final 7, 

• res,dir = resolved initial 7 and direct final 7, 

• res,frag = resolved initial 7 and fragmentation into final 7. 

Table 1: The cross section for the non-isolated (inclusive) photon, the isolated pho- 
ton, and the isolated photon -|- jet production. The photon transverse energy is in- 
tegrated over the range 5 < < 10 GeV. The initial electron and proton energies 
are E^. = 27.6 GcV and Ep = 920 GeV [49]. The results for the non-isolated photon 
are integrated over the whole range of y and r]'^. The isolated photon cross section 
is calculated for 0.2 < y < 0.7 and —1. < rf < 0.9. For the photon-|-jet production 
there are additional cuts, 4.5 GeV < E^^^ and -1. < 77^'^* < 2.3. 



[Pb] 


total 


dir,dir 


res,dir 


dir,frag 


res,frag 


Born 


0{as) 


box 


non-isolated 7 

E'^ cut 


240.6 


85.8 
(35.6%) 


5.0 
(2.1%) 


15.2 
(6.3%) 


60.0 
(25.0%) 


27.0 
(11.2%) 


47.6 
(19.8%) 


isolated 7 
y, Erp, 1]^ cuts 


37.77 


15.23 
(40.3%) 


1.76 

(4.7%) 


4.34 

(11.5%) 


12.63 

(33.4%) 


1.49 

(3.9%) 


2.33 
(6.2%) 


isolated 7+jet 
y, Ej,, 1]^ cuts 

Et\ Vjet cuts 


29.45 


11.60 
(39.4%) 


0.19 
(0.6%) 


3.41 
(11.6%) 


11.45 
(38.9%) 


1.20 
(4.1%) 


1.59 
(5.4%) 
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Table 2: The cross section for the non- isolated (inclusive) and isolated final photon 
with 5 < < 10 GeV. The initial energies are = 27.6 GeV and E.^ = 820 GeV [42]. 
The results for the isolated photon are obtained using the small cone approximation 
(Sec. 4) without and with cuts 0.2 < y < 0.9 and -0.7 < rj'^ < 0.9. 



[Pb] 


total 


dir,dir 


res,dir 


dir,frag 


res,frag 






Born 


0{as) 


box 








non-isolated 7 


222.0 


82.0 


4.8 


13.9 


54.8 


24.6 


42.0 


E^ cut 




(36.9%) 


(2.2%) 


(6.3%) 


(24.7%) 


(11.1%) 


(18.9%) 


isolated 7 


178.1 


82.0 


13.1 


13.9 


54.8 


5.1 


9.4 


Ej, cut 




(46.0%) 


(7.4%) 


(7.8%) 


(30.8%) 


(2.9%) 


(5.3%) 


isolated 7 


71.95 


23.60 


6.02 


6.53 


28.18 


2.34 


5.28 


y, Ej. cuts 




(32.8%) 


(8.4%) 


(9.1%) 


(39.2%) 


(3.3%) 


(7.3%) 


isolated 7 


35.34 


13.64 


3.26 


3.41 


11.88 


1.20 


1.93 


y, Ej,, Tj^ cuts 




(38.6%) 


(9.2%) 


(9.6%) 


(33.6%) 


(3.4%) 


(5.5%) 



Table 3: The cross section for the isolated final photon with E^ = 27.6 GeV, Ep = 
920 GeV, 0.2 <y < 0.9, -0.7 < 77^ < 0.9 and 5 < < 10 GeV [49]. The results 
are obtained with (upper row) and without (lower row) the small cone approximation 
(Sec. 5). 



[pb] 


total 


dir,dir 


res,dir 


dir,frag 


res,frag 


Born 


Oias) 


box 


approximated 


38.93 


15.23 


2.94 


4.34 


12.63 


1.49 


2.33 


exact 


37.77 


15.23 


1.76 


4.34 


12.63 


1.49 


2.33 
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Appendix A 

Kinematics and notation 



A.l 2^2 processes 

We study the production of photons and jets in the electron-proton scattering: 

e(Pe) piPp) e(pg) 7(p^) jet{pjet) X{px), (A.l) 

where four-momenta of particles are given in brackets, and assume that the mediating 
photon is quasi-real, 

Q'^-iPe-p'ef^O. (A.2) 

In such processes the emission of the mediating photon from the electron can be 
factored out (Sec. 2.3) and we can consider the photon-proton scattering: 

7(?) PiPp) l{P-i) jetiPjet) X{px), (A.3) 

where q — ypp is the four-momentum of the mediating photon. This photon may 
interact with a parton from the proton directly or as a resolved one, see Eq. (4.6). 
The corresponding 2^2 partonic processes are: 

a{Pa) b{Pb) c{Pc) d{pd), (A.4) 

where 

Pa = x^q = x^ype , Pb = xpp , Pc = Pj/z. (A.5) 
We introduce the standard variables: 

S=iPa+Pbf , t={pa-Pcf , U={pb-Pcf, (A.6) 

which are used in next Appendices and in formulae for the cross sections. 
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A.2. 2 -> 3 PROCESSES 



A. 2 2^3 processes 

Now, we consider the 2—^3 partonic processes contributing to the ep 67 (jet) X 
reaction: 7g — 'yqg (Fig. 2.3) and •yg •yqq (Fig. 2.4). We write 



and 



7(g) + q{p) -f{p^) + q{pi) + g{p2) 



7(g) + 9{p) 7(^7) + g(Pi) + g(P2 



(see Fie;. A.l) 



(A.7) 



We use the variables s, t and u defined in Appendix A.l, which for the direct both 



initial and final photon (x^ 



1) are given by 



u 



{p-p^f 



Finally, we define the scaled variables v and w. 

f = l + t/s , w = —u/{t + s). 



(A.9) 



(A.IO) 



The variables v and w are in the range from Otol: 0<f<l and < w < 1. 
Note, that for massless particles 

{pi+P2f = sv{l-w), (A.ll) 

and in the limit {pi +P2Y one obtains w ^ 1 {v can not be too low, since the final 
photon has a large transverse momentum). For the 2-^2 processes one has w = 1 by 
definition. 




7(^7) 



Figure A.l: An example of the 2 — > 3 NLO process contributing to the ep 
reaction. The four-momenta of particles are given in brackets. 



67 (jet) X 
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Appendix B 

Cross sections for 2^2 processes 



In this Appendix we have collected all the tree-level cross sections for 2—^2 processes 
which are included in our calculation. The virtual gluon corrections to the Born 
process are included in Appendix D together with the real gluon corrections. 

The cross sections for 2 ^ 2 processes shown in Figs. 2.1, 2.7-2.11 are given by 



the equation 

(5, t, u) = ■ |m"|2(s, t, u)-d{s + t + u), (B.l 

where [M^^^'^'^p are the squared matrix elements: 

\jr''-^'"'\\s,t,u) = -2(47r)Ve^(- + -), (B.2' 

^ s u 

|M"r(.,t,«) = -J(47r) W5(- + 7), (B.S: 

\W-'''\'is,t,u) = -^(47r) WJ(- + 7), (B.e; 

O S b 

I7ir™|2(s,i,^) = |7ir™|2(s,«,i), (B.7: 

|ir™|2(5,t,M) = {Anfaaseli- + -), (B.S; 

|M"r(.,t,x.) = (47r)^4 \-\{^ + '-) + ^1 , (B.g; 
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\M I [s,t,u) = \M I (s, 



■ai^ig 1 2 / 



\M I [s,t,u) = \M I [s,u,t), 



m-*Qg 1 2 / 



\M I (s,r,Mj = |M I (s, 



■gi^qg 1 2 / 



45^ + ^^ 

9 t2 



|M \[s,t,u) = \M \[s,u,t), 



|M"|2(s,t,«) = (47r)^a^ 



|M"|2(s,t,M) = (47r)2a1 



8 



4 , + 



9' ^2 



+ 



27 si 



4 ^2 + 

9^ 



\M' 



m^gg 1 2 , 



s, h) = (47r) dig 



1 , t -u, 
- + 7 



27^^ + 



3t2+M2 

8 s2 



8 ^2 + ^2 



t-u su st , 

^2 y2 ' 

The squared matrix element for the box process (Fig. 2.11) has a form: 

2 

2 



\M'''^^'\'{s,t,u) 



16\a J 



J2el] |Afn-.ii|' + |M22^22|'+ 



where the amphtudes M^j^ki (~ a"^) are given by Eq. (1) in Ref. [68], and 

Mi2^2l{s, t, U) = Mi2-^i2(s, U, t). 



2|Mn_22r + 2|Mi2_i2r + 2|Mi2_2ir) , (B.22) 
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Appendix C 
Three-body phase space 



The cross section for a 2 



3 process is given by a general formula: 



(C.l) 



where {PS)z is the three-body phase space in n dimensions: 



f d^pi r (i"p2 



(27r) 



n-l 



(27r)"-iy (27r)"-i 
• (2nr5^{q + p-p,-p,- p,)5-'{p'^)5-'(pI)5+(pI). 



(C.2) 



The integration of {PS)^ over all four-momenta of two final partons is discussed e.g. 
in [119, 18, 25]. However, in order to impose the isolation of the final photon as well 
as cuts for the jet, we need to restrict the momenta of the final partons. To achieve 
this goal we integrate {PS)^ considering various configurations of the final particles 
momenta, e.g. some collinear configurations, as described in Sec. 5.2. The integration 
of {PS)3 including a collinear configuration, where two final particles move (almost) 
parallel to each other, was previously discussed e.g. in [101]. In comparison with [101] 
we simplified the integration over the final energies. Below we derive all formulae for 
{PS)s which are applied in Appendix D to obtain the cross sections for the 0{as) 
corrections to the Born process including e.g. the configurations where a final parton 
is collinear with the final photon or with the initial electron or proton. 

The integration of {PS)^ as well as the calculations of the cross sections for the 
0{as) corrections to the Born process are performed in n — A — 2e dimensions using 
the dimensional regularization, see e.g. [1]. 

First, we integrate (C.2) over p2: 



(PSh = j^^s I dy.j dyiS^p'^)S^ipl)6^iiq + p-p,-p,n (C.3) 



Next, we choose the z-axis along the initial proton momentum in any frame of refer- 
ence, in which the collision is central, e.g. in the laboratory frame: 



p = xEp{l,...0,l), 



(C.4) 
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q = yE,{l,...Q,-l) (C.5) 

(the unspecified components are equal to zero). The colhsion has the rotational sym- 
metry and we can perform an integration over (n — 3) azimuthal angles of the photon. 
Performing this integration as well as the integration over the photon energy, and 
making the change of variables we obtain: 

/ dy^5+{p'^) = ^^^^^ JdvJ dw[vws(l - v)]-'sv. (C.6) 

Next, we chose the z and y axes this way that momenta p, q and are in the 2:7/- plane: 
p = x£'p(l, ...0, sina, coso;), (C.7) 

q — yEe{l, ...0, — sma, — cosa), (C.8) 

= E^{1, ...0,sma' ,cosa'), (C9) 

where a and a' are arbitrary, and the unspecified components in the additional {n — 4) 
dimensions are equal to zero. We write pi and p2 in this frame: 

Pi — E{1, ... cos ^3 sin ^2 sin ^1, cos ^2 sin ^1, cos ^1), (C.IO) 

P2^P + q-Pj-pi- (C.ll) 

Since the four-momenta p, q and p^ have only two non-zero space components, one 
can perform the integration over {n — 4) azimuthal angles of the final parton: 

J (TpiS+ipl) = ^^TZT) / / ^1 / ^2- (C.12) 

From (C.3), (C.6) and (C.12) one has 
d{PS)3 1 1 



dvdw 4(27r)"r(l - 2£) 



[vws{l — v)] ^sv- 



■J deisin^-^'eij de2sm-^'e2j dEE^-^'S+{{q+p-pj-pi)^). (CIS) 

Note, that the integration over E is straightforward as the argument of the S function 
is a linear function of pi: {q + p — P-y ^ PiY — {q + P ~ ^7)^ ~ 2pi(g + p — p^). 
The expression (C.13) was obtained with no approximations and can be used e.g. to 
calculate the cross section for 2 — > 3 processes in Part 5 of the phase space described 
in Sec. 5.2. 

In order to obtain {PS)s in a form suitable for considering coUinear configurations 
in Part 2, 3 and 4, we assume that 9i is small: 

9, ^ 0, (C.14) 
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APPENDIX C. THREE-BODY PHASE SPACE 



and take in Eqs. (C.7-C.9) a = tt or a = or a' = 0. With these assumptions we 
obtain: 

d{PS)s 



dvdw 



where 

• for a = tt: 

E^{1- w)yE,, (C.16) 

• for q; = 0: 

1 — vw 

• for a' — 0: 

E = "^^-"^ E,. (C.18) 

Of course, the same formulae for {PS)^ are vahd if we replace pi with 

P2 = E{1, ... cos 03 sin 02 sin 9i, cos6'2 sin 6'i, cos6'i), (C.19) 

Pl=P + q-P'y-P2- (C.20) 
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Appendix D 

Corrections of order 0{as) to the 
Born process 



In this Appendix we present formulae for 0{as) corrections (Figs. 2.2-2.4) to the 
Born process (Fig. 2.1). The general (unintegrated) formulae (Appendix D.l) and the 
formulae for the non-isolated photon production (Appendix D.2) are taken from the 
literature [18, 25]. The cross sections corresponding to various collinear configurations 
described in Sec. 5.2 are derived briefly in Appendices D.3, D.4 and D.5. 

D.l General formulae for 2^3 processes 

The general formula for the squared matrix element for the 75 ■jqg process has the 
form [18, 25]: 

010203616203 

+ 0{e% (D.l) 
where ft is an arbitrary mass scale, Cp — 4/3 and the terms Ai, Oj, 6j, c are given by 
Ai = aibi{al + bl) + 0262(02 + ^2) + ashial + 6^) (D.2) 

A2 = 0161(02 + 62) (03 + 63) - 0262(01 + 61) (03 + 63) 

-0363(01 + 61) (02 + 62) (D.3) 

A3 = albj + albl + ajbl (D.4) 

A4 = 01610262 + 01610363 -I- 0262O363 (D.5) 



and 
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APPENDIX D. CORRECTIONS OF ORDER 0{as) TO THE BORN PROCESS 



03 = PP2, (D.8) 
bi^Piq, (D.9) 
b2^PiPj, (D.IO) 
b3^PiP2^ -v{l-w), (D.ll) 

c^ppi. (D.12) 

The same formula (D.l) is valid for the squared matrix element of the '-/g '-/qq 
process, [aP'^^"''''' , but with the factor 1/2 instead of Cp, with no 0{e'^) term, and 
with the coefficients a[, a'2, a'^, b[, b'2, 63 and c' instead of ai, 02, 03, 61, 62, ^3 and c, 
where the coefficients with the prime are given by: 

a'l = -P2q: (D.13) 

4 = -P2P7, (D.14) 
4 = a3, (D.15) 
b[ = b,, (D.16) 
6^ = 62, (D.17) 
6'3 = -c, (D.18) 

c' = -63. (D.19) 

The above squared matrix elements are used (with £ = 0) to calculate the cross sections 
for the 0{as) corrections to the Born process in Part 5 of the phase space (see Sec. 
5.2): 

d'^p-y svTT 2s dvdw 
where the three-body phase space is given by Eq. C.13. 
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D.2. INCLUSIVE PHOTON CROSS SECTION 



D.2 Inclusive photon cross section 

The cross sections for the 0{as) corrections to the Born process, integrated over all 
momenta of final partons are given by [18, 25]: 



9(1 -w 



2 q 



ci6{l-w) + 



C2 



1 — w 



e3 ,5141^1 +fefl,i(l-«,) + 



1 — w 



+ 



cW^l—^ — ) + cAlia-^ + c^lnv + celnil - vw)+ 
.l-wj+ J 11'^ 



C7 ln(l — V + vw) + Cg ln(l — + Cg In w + Cio ln(l — w) + Cii+ 



\n.{l — V -\- vw) Inw ln((l — i>w)/(l — ^;)) 
C12 ; 1- C13- I-C14 



1 — w 



1 — w 



1 — w 



(D.21) 



= 6(1 - w) -e 



a^as 4I 



ITS 



2 1' 



C4 In — + C5 In V + Cg ln(l — vw)+ 



1^ 



C7ln(l — V + vw) + C8ln(l — v) + c'glnw + c'ioln(l — w) + c'n + 



I ln(l — V + vw) , 

^12 ^ Z I" "^IS" 



1 



w 



Inw , ln((l — t'w)/(l — v)) 
; I" ^14 ; 



(D.22) 



where the coefficients c and c' are given in Appendix B in the paper [18]. We use 
formulae (D.21) and (D.22) to obtain the 0{ols) corrections for the non-isolated photon 
production (Chapters 3, 4) as well as to obtain the Oiots) corrections for the isolated 
photon in Part 1 of the phase space, as defined in Sec. 5.2 (Chapters 5, 6). 



D.3 Collinear configuration (^i||^e or V'2\\V(!) 

In this Appendix we consider the configuration corresponding to the Part 2 of the 
phase space defined in Sec. 5.2. First, we assume that the momentum of the final 
quark, pi, is almost parallel to the momentum of the initial photon originating from 
the electron, q {— ype). We orient the axes this way that 

q = yE,{l,..., 0,0,1), (D.23) 
p^xEpil,..., 0,0,-1), (D.24) 
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APPENDIX D. CORRECTIONS OF ORDER 0{as) TO THE BORN PROCESS 



p-y = E^{1, 0, sin a' , cos a'), (D.25) 

Pi = E{1, cos6'2 sin 9i, cos 6*1), (D.26) 

P2 = q + p-p^-pi, (D.27) 
where 

= xEp{l -v)+ yEeVw, (D.28) 



— v)vw 

sin a' = ^ \^ — , (D.29) 
E^ 

, xEJl — v) — vEeVW 

cos a' = ^ ^ — — — . (D.30) 

Ej 

E = {l-w)yE,, (D.31) 

< ^1 < e,ut < 1. (D.32) 

From Eqs. (C.l, C.15, C.16, D.1-D.19) we obtain the cross sections which contain 
terms ~ 1/e being singular in 4 dimensions: 

E^^SS1\ = J-.l^i^|M^™|2~0(i) + 0(l) + ..., (D.33) 

d'^P-y \ singular SVTT 2s dvdw £ 



d^P^ \ singular SVTT 2s dvdw S 

We remove these singularities applying the standard factorization procedure (see e.g. 
[18, 25]): 

E,^g^ = ^7^^ I - / diH,,{i,^j:)E,'^^ {is^M, (D.35) 

a-^P^ d'^P^ \ singular J O'^P-y 

^7^^ = ^7^^ I - / d^H,,{^,^^)E,— {^s,^t,u), (D.36) 

a'^P-y d-^Pj \ singular i d^^p^ 

with 
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D.4. COLLINEAR CONFIGURATION {Pi\\Pp OR P2||i3^) 



-C/-V — {s,t,u) — 



3s 



r(l -e) \ tu 



, 't u 
^ - + 7 



-2e 



;d.38) 



s, t, u) — 



nfi s 



3s 2 r(i - £) V 



.s t 



2e 



;d.39) 



where ft is an arbitrary mass scale. Finally, we obtain the singular-free 0{as) cor- 
rections to the Born process in the region of the phase space labeled as Part 2 (Sec. 
5.2): 



da. 



d^p. 



ns^vw v{l — v) 



;d.4o) 



daZ' 



"S(2) 



= ^(1 - w) 



a^ase^q 1 1 



[w^ + (1 - In 



+ 1 



(D.41) 



These cross sections (D.40, D.41) are equivalent to the results quoted in Sec. 5.2 in 
Eqs. (5.10, 5.11). 

The squared matrix element iM^^^^'^^p calculated for the gluon moving (almost) 
parallel to the initial electron, P2I I? = yPe (see the notation in Eq. (A. 7)), contains no 
coUinear singularities and the cross section da'^^^^'^^ corresponding to this configura- 
tion can be obtained performing the exact numerical calculations as in Part 5 of the 
phase space (see Appendix D.l) or it can be neglected as it is a term of order 0{9^ 



cut) 



On the other hand, the configuration corresponding to the final antiquark coUinear 
to the initial electron, p2\\q = yPe (see the notation in Eq. (A.8)), leads to the 
same contribution in the cross section da^jy^'^'^ as the quark coUinear to the electron, 
Pill? — yPe: a-nd this is included in Eq. (5.5) in the summation over 2Nf fiavours. 



D.4 CoUinear configuration (^i||^ or P2||Pp) 

Herein we derive the cross sections da'l''^^'^^ and da^J'^^'^'^ corresponding to the region 
of the phase space labeled in Sec. 5.2 as Part 3. First, we consider the final gluon 
moving almost parallel to the initial quark and the final antiquark moving almost 
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parallel to the initial gluon, P2\\p = xpp (see the notation in Eqs. (A. 7, A. 8)). We 
write the four-momenta in the laboratory frame with the z axis chosen along the 
direction of the initial proton momentum: 

q = yE,{l,..., 0,0,-1), (D.42) 

p = xEp{l,... ,0,0,1), (D.43) 

— E^{1, ...,0, sin a', cos a'), (D.44) 

P2 = E{1, ...,cos^2sin^i,cos^i), (D.45) 

Pi^q + P-Pj-P2, (D.46) 
where 

E^ = xEp{l -v)+ yE^vw, (D.47) 



Jsil — v)vw 

sin a' = " \. , (D-48) 

, xEJl — v) — yEgVW ,^ 

cos a' = — ^ — — — , (D.49) 

E^ 

E =. '-^^xEp, (D.50) 
1 — vw 

0<9i< 9,ut < 1- (D.51) 

Prom Eqs. (C.l, C.15, C.17, D.1-D.5, D.13-D.19) we have the cross sections which are 
singular in 4 dimensions: 

E^^SSll = J-.l4^|]r™|2^0(l) + O(l) + ..., (D.52) 

d'^p-y \ singular SVTT 2s dvdw S 

E,'^-^SS1\ = J-.l^^|iir™r^O(i) + 0(l) + ..., (D.53) 

d'^P^ \ singular SVTT 2s dvdW 8 

As in Appendix D.3, we remove these singularities applying the factorization proce- 
dure: 

^7^^ = ^7^^ I - / diH,,{i,^j)E,^^ {is,t,iu), (D.54) 

d'^p-y d-^p.y \ singular J 0"^^^, 
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da 



-C/-V — -C/- 



d'^Py d^Py \ singular 



with 



as 
Gaet 



(D.56) 



and 



H(^n)- 1 + ^' l^^vrAn^ r(l^.) 



s r(i -£) \ tu 



(D.57) 



-(l-5)(- + ^)-2s 
s u 



;d.58) 



The final singular-free formulae for the 0{as) corrections to the Born process in Part 
3 of the phase space are: 



2 4 

7rs^(l — v) 



Cp ( WW + 



1 



V 



vwj 



1 — X 



(D.59) 



da 



as (3) 



= B 1 — — viy H 



+ (l-x 



where 



1 



a; = 



1 — i'w 



(D.60) 



(D.61) 



These formulae are equivalent to the formulae (5.12, 5.13) in Sec. 5.3. 

The squared matrix element [TW^'^^^'^^p calculated for the final quark moving par- 
allel to the initial quark, pi\\p = xpp (see Eq. (A. 7)), contains no coUinear singularities 
and the cross section da^^^^'^^ corresponding to this configuration can be either cal- 
culated numerically or neglected. 

The configurations corresponding to the final quark or antiquark coUinear to the 
initial electron, pi\\p = xpp or P2\\p = xpp (see Eq. (A.8)), give the same contributions 
to the cross section daZ^J^^V'^ and this is included in Eq. (5.5) in the summation over 



2A^y^ flavours. 



'as (3) 
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D.5 Collinear configuration {pi\\p^) 

Now, we choose the z axis in the direction of the final photon momentum: 

q = yEg{l, 0, — sina, — cosa), (D.62) 

p — xEp{l, ...,0,sma,cosa), (D.63) 

= 0,0,1), (D.64) 

Pi — E{1, ...,cos^2sin^i,cos^i), (D.65) 

P2 = g + P -P7 -Pi, (D.66) 

with 

E^^xEp{l-v)+yEeVW, (D.67) 



sin a = ^ ' , (D.68) 

E^ 

cosa = ^^2^4^^^, (D.69) 
E^ 

E^ "^^-^^ E, ^-'-±^E,. (D.70) 

l-v + vw ^ t + u ^ ^ ^ 

We assume that the angle between the final quark momentum and the final photon 
momentum is small: 

< ^1 < 5e„t < 1, (D.71) 

where Scut is a cut-off parameter. The partonic cross sections derived from Eqs. (C.l, 
C.15, CIS, D.1-D.12, D.20) contain the 1/e singularity: 

E,^-^^\ = J-.l^i^|ir™r~0(i) + 0(l) + ..., (D.72) 

d'^Pj \ singular SVTT 2s dvdw S 

^^^Pr, = J_.l^i^|Ar™|^~0(l) + O(l) + ... (D.73) 

a'*P-y \singular SVTT 2s dvdw S 

The singularities are removed in the factorization procedure: 

^ "^a5(4) ^ 

d'^p^ 
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singula 



dCH,,{C, i,)E,^^{s, t/C, u/0, (D.74) 



as (4) 



= E. 



singular 



r da'^^^^^ 
J d^H^g{tfi)E^—j^{s,t/tu/0, 



d^p^ 



with 



la 2 1 + (1 -OV^T^A^V r(i 



J r(i-2£) 



+ 0{e), 



-B-v — {s,t,u) = 8E^—^ — (•5,^,1*), 



^ d^p^ 



£/-v 77^ [s,t,u) — -E^ — [s,t,u), 



d^p^ 



8 ^ d^p^ 



(D.75) 

(D.76) 
(D.77) 
(D.78) 



where the partonic cross sections da'^'^^'^^ and da^'^^'^^ are given in Appendix D.4. 
Finally, the 0{as) corrections to the Born process in Part 4 of the phase space are 



E^ 



daZ" 



^■yqg 



as (4) 



d^p^ 



a^ase'*^ 1 + (1 _ -0) 



1-v 



In 



+ z 



(D.79) 



d^p-. 



a^asefj 1 v"^ + (1 - vf 



TTS'^Z 2 ^(l — v) 



l + {l-z)\ E'JUl-z) 



In 



+ z 



where 



z — 1 — V -\- vw, 
vw 

V 



(D.80) 

(D.81) 
(D.82) 



(D.83) 



1 — V + vw 
Since 5 cut is small, we have 

where Rcut is the cut-off radius in the rapidity and azimuthal angle space, see Sec. 5.2. 
From Eqs. (D.79), (D.80) and (D.83) one obtains Eqs. (5.15, 5.17) corresponding to 
Part 4 of the phase space as well as Eqs. (4.20, 4.22) for the subtraction terms in the 
small cone approximation. 
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D.6 Collinear configuration (P2IIP7) 

Finally, we assume that the momentum p2 is almost parallel with the momentum of 
the final photon, p^: 

q — yEg{l, 0, — sin a, — cos a), (D.84) 

p = a:£'p(l, 0, sina, cosa), (D.85) 

= 0,0,1), (D.86) 

P2 = E(l, ...,cos6'2sin6'i,cos6li), (D.87) 

Pi^q+P-Pj-P2- (D.88) 

The calculation of the cross section for the process — > ^qq is the same as in the 
previous Appendix D.5, and the formula for da'^J^^'^'^ is given by Eq. (D.80). This 
is included in the cross section for the 7p collision in the summation over quarks and 
antiquarks {2Nf flavours), see Eq. 5.5. 

The calculation for the process 75 — > ''^qg with the gluon momentum (^2) parallel 
to p^ is different than the corresponding calculation in Appendix D.5 for the quark 
momentum (pi) parallel to p-y, since now we deal with the cross section which contains 
no singularities. 

The cross section da'^J^^'^^ is equal to 

d'^p^ svn 2s dvdw 

where {PS)^ and [M^'' "''^^\'^ are given in Eqs. (C.15) and (D.l), respectively. This 
cross section contains no singularities and is proportional to O^^^. In the Part 4 of the 
phase space we assume that w < Wcut- This cut-off restricts the phase space of the 
final gluon this way that the gluon can not be too soft, so the cross section (D.89) 
contains no large terms due to the emission of soft gluons and is small ~C^(^cwt)- 
can be either calculated numerically or neglected if 9cut is sufficiently small. 

Similar configuration (gluon moving parallel to the final photon) was considered 
in Chapter 4, when the subtraction cross section in the small cone approximation 
was discussed. In this case we had no cut-off for w, and large terms due to the soft 
gluon emission might contribute to the cross section, so the subtraction term da]^^^ 



(4.21) could not be neglected. We have obtained the formula (4.21) for the subtraction 
term from Eqs. (C.15, D.l, D.89) keeping all terms of order C(6'^„J and neglecting 
higher powers of 9cut- 
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